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FOREWORD 


This  document  is  Volume  XX  of  the  Interim  Report  series  for  the  Passive  Nosetip  Tech- 
nology (PANT)  program.  A summary  of  the  documents  in  this  series  prepared  to  date  is  as  fol- 
lows : 


Volume  1 
Volume  II 
Volume  III 


Volume  IV 


Volume  V 


Volume  VI 
Volume  VII 

Volume  VIII 

Volume  IX 


- Program  Overview  (U) 

- Environment  and  Material  Response  Procedures  for  Nosetip  Design  (U) 

- Surface  Roughness  Data 

Part  I — Experimental  Data 

Part  II  - Roughness  Auginented  Heating  Data  Correlation  and  Analysis  (U) 
Part  III  - Boundary  Layer  Transition  Data  Correlation  and  Analysis  (U) 

- Heat  Transfer  and  Pressure  Distributions  on  Ablated  Shapes 
Part  I - Experimental  Data 

Part  II  - Data  Correlation 

- Definition  of  Shape  Change  Phenomenology  from  Low  Temperature  Ablator 
Experiments 

Part  I —Experimental  Data,  Series  C (Preliminary  Test  Series) 

Part  II  —Experimental  Data,  Series  D (Final  Test  Series) 

Part  III  - Shape  Change  Data  Correlation  and  Analysis 

- Graphite  Ablation  Data  Correlation  and  Analysis  (U) 

- Computer  User's  Manual,  Steady-State  Analysis  of  Ablating  Nosetips 
(SAANT)  Program 

- Computer  User's  Manual,  Passive  Graphite  Ablating  Nosetip  (PAGAN)  Pro- 
gram 

- Unsteady  Flow  on  Ablated  Nosetip  Shapes  - PANT  Series  G Test  and  Analy- 
sis Report 


Preceding  page  blank 


i ii 


Volume 

Volume 

Volume 

Volume 

Volume 

Volume 

Volume 

Volume 

Volume 

Volume 

Volume 

Volume 

Volume 

Volume 


X - Summary  of  Experimental  and  Analytical  Results 

XI  —Analysis  and  Review  of  the  ABRES  Combustion  Test  Facility  for  High  Pres- 

sure Hyperthermal  Reentry  Nosetip  Systems  Tests 

XII  - Nosetip  Transition  and  Shape  Change  Tests  in  the  AFFDL  50  MW  RENT  Arc  - 

Data  Report 

XIII  - An  Experimental  Study  to  Evaluate  Heat  Transfer  Rates  to  Scalloped  Sur- 

faces - Data  Report 

XIV  - An  Experimental  Study  to  Evaluate  the  Irregular  Nosetip  Shape  Regime  - 

Data  Report 

XV  — Roughness  Induced  Transition  Experiments  - Data  Report 

XVI  — Investigation  of  Erosion  Mechanics  on  Reentry  Materials  (U) 

XVII  —Computer  User's  Manual,  Erosion  Shape  (EROS)  Computer  Program 

XVI II  - Nosetip  Analyses  Using  the  EROS  Computer  Program 

XIX  — Hydrometeor/Shock  Layer  Interaction  Study 

XX  — Investigation  of  Flow  Phenomena  Over  Reentry  Vehicle  Nosetips 

XXI  - Flight  Implications  of  Low  Temperature  Ablator  Shape  Data  (U) 

XXII  - Coupled  Erosion/Ablation  of  Reentry  Materials 

XXIII  - Reentry  Vehicle  Nosetip  Response  Analyses 


This  report  was  prepared  by  Aerotherm  Division/Acurex  Corporation  under  Contract  F04701- 
71 -C-00E7.  Volumes  I through  IX  covered  PANT  activities  from  April  1 9/4  through  April  1975. 
Volumes  X through  XV  represent  contract  efforts  from  May  1973  to  December  1974.  Volumes  XVI 
through  XVIII  describe  the  background,  development,  and  check  out  of  the  PANT  EROsion  Shape 
(EROS)  computer  code.  These  volumes  document  efforts  performed  under  supplementary  agreements 
to  the  Minuteman  Natural  Hazards  Assessment  program  (Contract  F04701 -74-C-0069 ) between  April 
1974  and  March  1975.  Volumes  XIX  through  XXIII  document  additional  analyses  performed  between 
December  1974  and  June  1975. 


This  work  was  administered  under  the  direction  of  the  Space  and  Missile  Systems  Organiza- 
tion with  Lieutenant  A.  T.  Hopkins  and  Lieutenant  E.  G.  Taylor  as  Project  Officers  with  Mr.  W. 
Portenier  and  Dr.  R.  L.  Baker  of  the  Aerospace  Corporation  serving  as  principal  technical 


monitors. 
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Project  Officer 
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SECTION  1 


INTRODUCTION 

During  the  period  April  1971  through  June  1975,  the  Aerotherm  Division  of  Aairex  Corpora- 
tion was  involved  ir  the  development  of  reentry-vehicle,  passive  nosetip  design  tools  under  the 
SAMSO  sponsored  Passive  Nosetip  Technology  (PANT)  program.  The  effort  included  experimental 
and  analytical  investigations  of  the  fundamental  phenomena  controlling  the  response  of  nosetip 
materials.  Results  of  the  first  phase  of  this  program  were  reported  in  the  PANT  program.  In- 
terim Final  Report,  Volumes  I through  VIII,  dated  January  1974  (Reference  1-1).  The  signifi- 
cant findings  consisted  of  the  following: 

• The  importance  of  the  nosetip  surface  roughness  on  boundary  layer  transition  ana 
heat  transfer  was  demonstrated  through  wind  tunnel  simulation  experiments. 

• Wind  tunnel  simulation  experiments  also  demonstrated  that,  the  surface  temperature 
to  boundary  layer  edge  temperature  ratio  significantly  Influences  the  location  and 
occurrence  of  boundary  layer  transition. 

• Nosetip  ablation  tests  in  both  wind  tunnel  and  hyper'.hermal  facilities  indicated 
that  surface  roughness  patterns  (i.e.,  scallops)  which  affect  heat  transfer  develop 
during  turbulent  ablation. 

• Consideration  of  roughness  effects  on  heat  transfer  indicated  that  the  micro:, iechani  - 
cal  particulate  component  of  graphitic  nosetip  mass  loss  is  negligible  and  that  gra- 
phite mass  loss  and  surface  temperature  data  can  be  predicted  utilizing  equilibrium 
thermochemistry  theory. 

• Nosetip  shape-change  data,  primarily  from  wind  tunnel  tests,  indicated  that  nosetip 
ablation  can  produce  irregular  (slender  or  asynmetric)  shapes  which  are  prone  to 
gross  failures. 

The  primary  objective  of  the  second  phase  of  the  PANT  program  was  to  define  the  reentry 
regimes  which  produce  irregular  shapes.  This  phase  covered  the  period  from  May  1973  through 
December  1974  and  was  comprised  of  both  analytical  and  experimental  efforts.  Results  of 


the  second  phase  are  documented  in  the  Interim  Report,  Volumes  IX  through  XV.  The  most  impor- 
tant results  from  those  efforts  are: 

• Shape  change  response  of  flight  materials  in  the  AFFDl  5G  MW  arc  jet  was  consistent 
with  the  response  of  numerous  tests  of  low  temperature  ablator  models  tested  in  the 
NSWC  wind  tunnel . 

• Scallops,  which  form  on  both  LTA  models  in  wind  tunnels  and  graphite  models  in  the 
50  MW  arc  jet,  cause  the  convective  heat  transfer  rate  to  increase  over  the  corre- 
sponding smooth  surface  predictions.  Flowever,  the  roughness  augmentation  factor 
iratio  of  actual  to  smooth  wall  heat  transfer  rate)  does  not  scale  with  free  stream 
Reynolds  number,  in  contradistinction  to  the  behavior  exhibited  by  surfaces  having 
sand  grain  type  roughness. 

• The  Reynolds  number  condition  for  the  development  of  irregular  shapes  in  wind  tunnel 
low  tempei  ature  ablator  tests  varies  significantly  with  stream  total  temperature, 
initial  nose  radius,  and  surface  roughness.  During  the  wind  tunnel  tests,  the  for- 
mation of  ■frregular  shapes  was  accompanied  by  a significant  increase  in  model  vibra- 
tional accelerations  associated  with  unsteady  or  quasi -unsteady  flow  phenomena. 

• The  transition  criteria  developed  in  the  first  phase  of  the  PANT  program  were  extended 
to  flight  conditions  by  including  the  dependence  on  ablation  and  species  dissociation. 
Wind  tunnel  tests  demonstrated  that  the  criteria  are  applicable  over  a range  of  body 
sizes  and  shapes.  Allowing  for  material  surface  roughness  uncertainty,  the  extended 
criteria  agree  with  available  flight  nosetip  transition  data. 

<•  Shape  change  predictions  were  compared  with  flight  data  for  seven  flights.  Modifi- 
cations to  the  modeling  and  solution  procedures  improved  the  predictive  capability, 
but  the  prediction  of  the  irregular  shape  regime  remained  uncertain. 

One  of  the  most  important  results  of  the  second  phase  was  the  determination  of  the 
importance  of  the  irregular  shape  regime  to  nosetip  performance.  In  addition,  the  sensitivity 
of  shape  change  response  predictions  to  modeling  of  the  aerothermal  environment  was  identified 
as  an  important  factor  in  limiting  the  current  capability  to  predict  the  onset  of  irregular 
shapes  in  flight  predictions. 

The  third  phase  of  the  PANT  program  was  directed  toward  better  characterizing  the  depen- 
dence of  Irregular  shape  development  on  the  aerothermal  environment  and  identifying  needed 
improvements  in  analytical  modeling  to  better  predict  nosetip  shape  response  in  flight 
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environments.*  Specific  areas  which  were  addressed  included  coupled  ablation/erosion  effects, 
inviscld  and  boundary  layer  flow  effects,  nosetip  transition  on  advanced,  "smooth"  materials, 
and  irregular  shape  regime  onset.  This  report  documents  the  results  oi  efforts  on  the  inviscid 
flow,  "smooth"  wall  nosetip  transition,  and  turbulent  boundary  layer  heat  transfer  tasks. 

These  three  tasks  are  documented  in  Sections  2,  3,  and  4,  respectively,  each  of  which  is  self- 
contained.  The  most  important  conclusions  and  recommendations  resulting  from  these  tasks  are 
summarized  in  Section  5. 

These  three  tasks  are  documented  together  because  of  their  inherent  interdependence; 
the  subject  phenomena  are  the  base  of  the  environment  which  the  nosetip  experiences.  Particu- 
lar areas  of  concern  at  the  initiation  of  this  phase  were; 

• Inviscid  Flow  - The  shock  shape  prediction  is  important  in  evaluating  entropy  layer 
effects  on  turbulent  heat  transfer  rates  and  in  evaluating  the  effect  of  hydrometeor 
impact  on  nosetip  material  erosion.  Numerical  sensitivity  studies  indicated  that 
the  existing  shock  shape  prediction  procedure  is  inaccurate  for  certain  shapes,  and 
and  the  objective  of  the  inviscid  flow  field  task  was  to  evaluate  the  feasibility 
of  implementing  a procedure  based  on  solving  a direct,  thin  shock  layer  model  in 
order  to  improve  the  shock  shape  predictions. 

« Transition  - In  order  to  provide  a technique  for  Improved  prediction  of  transition 
on  finer  micros tructu re  graphite  nosetips,  the  earlier  transition  studies  were 
extended  to  generalize  the  PANT  rough  wall  transition  criteria  to  account  for  the 
effect  of  free  stream  disturbances  and  to  begin  characterization  of  the  free  stream 
disturbances  to  which  nosetips  are  subjected  in  the  atmosphere. 

» Turbulent  Heat  Transfer  - In  the  shape  change  codes,  smooth  wall  turbulent  heat 
transfer  prediction  is  based  on  solutions  of  the  integral  momentum  equation  for 
the  momentum  thickness,  a correlation  of  skin  friction  as  a function  of  momentum 
thickness  edge  Reynolds  number,  and  application  of  Reynolds  analogy  to  obtain  the 
heat  transfer  coefficient  from  the  skin  friction.  Nonsimilar  effects  resulting 
from  streamwise  pressure  gradient  and  boundary  layer  edge  normal  entropy  gradients 
are  neglected,  Surface  microroughness  effects  are  accounted  for  with  an  empirical 
multiplicative  augmentation  factor,  and  macroroughness  (scallop)  effects  are  eval- 
uated using  the  same  empirical  correlation,  wish  the  value  of  effective  surface 
roughness  determined  by  comparing  predictions  with  ground  test  and  flight  data. 
Comparisons  with  ground  test  data,  particularly  ballistic  range  data,  indicated 
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thit  there  are  inadequacies  in  the  procedure  which  may  be  quite  important  in  flight 
predictions.  The  objective  of  this  task  was  to  identify  areas  of  inadequacy;  par- 
ticularly with  respect  to  streamwise  pressure  gradient  and  normal  entropy  gradient 
effects,  and  to  modify  the  procedure  to  eliminate  the  inadequacies. 

SECTION  1 REFERENCES 

1-1.  "Passive  Nosetip  Technology  'PANT)  Program,  Interim  Report,  Volumes  I through  VIII," 

SAMSO-TR-74-86,  Aerotherm  Report  74-90,  Aerotherm  Division,  Acurex  Corporation,  January 
1974. 
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INVISCIO  FLOW  MODELING  STUDIES 
by 

Larry  Cooper 
Peter  Overly 


2-1 


list  of  symbols  for  SECTION  2 


A 

9 

H 

o 

m 

M 

An 

AN 

P 

Pt 

rs 

Ar 

7 

R 


R* 


T 


t 


w 

xs 

a 

6 

A 


Y 


area 

gravitational  constant 

stagnation  point  shock  stand-off  parameter 

mass  flow  function  (see  Equation  (2-10)) 

Mach  number 

flow  increment  spacing  along  nonrial 

local  shock  layer  thickness  measured  perpendicular  to  the  body 

static  pressure 

total  pressure 

radial  coordinate  of  shock 

radial  flow  increment  in  free  stream 

mean  radial  distance 

radius  of  curvature 

gas  constant 

total  temperature 

mass  flow  rate 

axial  coordinate  of  shock 

difference  between  body  angle  and  local  flow  deflection  angle 

flow  deflection  angle 

stagnation  point  shock  stand-off  distance 

ratio  of  specific  heats 
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LIST  OF  SYMBOLS  FOR  SECTION  2 (Concluded) 


p density 

9g  body  angle 

9$  shock  angle 

Subscripts 
B body 

i update  index  on  flow  angularity 

j index  along  normal  starting  at  body 

k index  of  normal  location  along  body 

NF  new  flow 

st  evaluated  at  stagnation  point 

s evaluated  at  shock 

w evaluated  at  wall 

1 frees tream 


2 


downstream  of  shock 


2.1  INTRODUCTION  AND  OBJECTIVES 


Reentry  of  a strategic  missile  into  the  earth's  atmosphere  at  high  velocity  produces 
evere  heating  loads  on  the  vehicle  nosetip.  As  the  vehicle  passes  through  the  atmosphere 
v.e  air  is  processed  by  means  of  a shock  wave  formed  ahead  of  and  around  the  body.  This 
Shirk  wave  raises  the  temperatures  and  pressures  adjacent  to  the  body  surface.  Meanwhile, 
next  to  the  body  surface,  viscous  effects  create  a dissipation  of  kinetic  energy  which  is 
converted  into  heat  ana  is  convected  to  the  surface.  The  driving  potential  for  this  convec- 
tive heat  transfer  through  the  boundary  layer  is  proportional  to  the  mass  flow  per  unit  area 
at  the  \.<ne  of  the  boundary  layer  and  the  convective  heat  transfer  coefficient.  The  mass 
flow  per  rit  area  is  a direct  function  of  the  local  edge  conditions.  The  total  thermodyna- 
mic temperature  and  pressure  at  the  boundary  layer  edge  are  determined  by  the  local  inclina- 
tion of  the  .blique  shock  through  which  that  particular  streamline  passed.  In  addition,  the 
convective  heav  transfer  coefficient  can  be  related  to  the  local  thickness  of  the  boundary 
layer.  This  boundary  layer  growth  is  affected  by  the  local  edge  venditions  all  along  the 
boundary  layer  edge.  Hence,  the  convective  heat  transfer  coefficient  can  be  related  to  the 
integrated  effect  of  boundary  layer  edge  conditions.  Put  another  way  the  convective  heat 
transfer  coefficient  is  related  to  the  entire  shape  of  the  shock  from  the  stagnation  stream- 
line to  the  point  where  the  local  edge  streamline  crosses  the  shock.  Hence,  it  is  apparent 
that  both  the  shock  shape  and  local  slopes  need  to  be  defined  accurately  in  order  to  perform 
meaningful  heat  transfer  and  ablation  calculations. 

The  nosetip  of  the  vehicle  assumes  a myriad  of  shapes  during  the  trajectory.  Each  in- 
stantaneous shape  of  the  nosetip  and  heat  shield  produces  its  own  shock  shape.  This  shock 
shape  produces  a unique  set  of  edge  conditions  for  driving  the  heat  transfer  and  ablation 
phenomena,  thereby  producing  a unique  shape  change.  Thus,  the  entire  shape  change  history 
is  dependent  on  the  shock  shape,  and  vice  versa.  It  is,  therefore,  required  that  a means  of 
accurately  predicting  shock  shape  for  arbitrarily  shaped  bodies  be  inherent  to  a good  shape 
change  analysis. 

The  means  by  which  the  shock  shape  and  boundary  layer  calculations  are  related  is  the 
so-called  entropy  swallowing  calculation  which  is  performed  to  determine  the  local  boundary 
layer  edge  conditions.  This  is  accomplished  by  matching  the  mass  flow  in  the  boundary  layer 
to  an  equal  amount  of  flow  passing  through  the  bow  shock.  Thus,  by  flow  rate  balance  at 
each  point  on  the  surface,  one  can  determine  the  point  at  which  the  local  boundary  layer  edge 
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streamline  pa-sed  through  the  bow  shock.  Knowing  the  local  inclination  of  the  bow  shock 
thus  determines  the  total  conditions  on  the  local  edge  streamline.  Given  the  local  static 
pressure  on  the  body,  the  local  boundary  layer  edge  conditions  can  be  determined. 

A second  importance  can  be  attached  to  the  determination  of  slock  shape  with  regard 
to  possible  reentry  flight  through  weather  or  dust  clouds.  The  presence  of  the  shock  wave 
ahead  of  the  nosetip  acts  as  a means  of  reducing  the  destructive  effects  of  particle  impact 
due  to  particle  demise  and  particle  deflection  behind  the  shock.  From  this  standpoint,  accu- 
rate knowledge  of  the  shock  standoff  all  around  the  nosetip  is  important.  In  addition,  if 
the  shock  shape  method  also  gives  information  about  the  flow  field  between  the  shock  and  the 
body,  that  is,  in  the  "shock  layer",  then  flow  deflection  and  part  cle  demise  calculations 
can  be  performed. 

Previous  means  of  determining  the  shock  shape  within  the  PANT  shape  change  codes  de- 
pended upon  correlative  or  semicorrelative  techniques.  The  earliest  shock  shape  determina- 
tion method  was  a correlation  of  shock  shape  versus  local  body  angle  for  spheres  in  hyper- 
sonic flow  (Reference  2-1 ) . Obviously,  such  a procedure  is  not  applicable  to  a wide  variety 
of  body  shapes  and  flow  conditions.  Later,  the  method  was  improved  (Reference  2-2)  by  incor- 
porating a local  shock  expansion  technique  In  the  supersonic  region.  However,  in  the  sub- 
sonic regions  the  older  technique  based  upon  local  body  angle  was  used.  In  its  overall 
ability  to  determine  shock  shape,  this  previous  method  was  shown  to  be  a great  improvement 
over  the  older  method.  Realistically,  the  improved  method  still  had  three  basic  weaknesses: 
(1)  for  blunted  shapes  that  have  large  subsonic  regions  the  body  angle  correlation  is  inade- 
quate; (2)  the  stagnation  point  shock  standoff  must  be  determined  independently,  and  (3)  the 
applicability  of  the  general  shock  expansion  method  for  a large  class  of  arbitrary  shapes  is 
questionable.  That  Is,  the  method  lacks  a certain  amount  of  generality  In  its  ability  to 
cope  with  some  classes  of  shapes  encountered  in  nosetip  applications. 

The  method  to  be  described  in  this  report  is  intended  to  be  free  of  the  shortcomings 
of  the  previous  shock  shape  prediction  methods  discussed.  It  is  applicable  to  both  the 
subsonic  and  supersonic  regions.  Includes  determination  of  the  stagnation  point  shock  stand- 
off, and  is  completely  general  in  its  applicability  to  arbitrary  shapes.  In  addition,  the 
method  predicts  the  approximate  shock  layer  flow  field  if  such  information  is  required  for 
other  purposes. 
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ANALYTICAL  APPROACH 


2.2.1  Thin  Shock  Layer  Theory 

Exact  numerical  solutions  to  the  problem  of  blunt  body  shock  layer  flow  are  time  con- 
suming and  expensive  from  a computational  standpoint.  For  hypersonic  flows,  > 5,  the  be- 
ht  uf  the  flow  field  is  such  that  more  approximate  and  less  expensive  techniques  can  be 
employed  to  study  shock  shapes  about  axisymmetric  bodies.  This  approach  is  known  as  thin 
shock  layer  theory  (References  2-3,  2-4,  2-5,  2-6,  and  2-7).  The  essence  of  the  method  lies 
in  the  simplifications  that  can  be  made  when  it  can  be  assumed  that,  generally,  the  flow  in 
the  shock  layer  is  "along"  the  body,  as  illustrated  in  Figure  2-1.  What  is  meant  here  is 
that  the  flow  is  nearly  perpendicular  to  local  body  normals  to  the  surface  and  that  the  flow 
curvati re  is  in  the  same  sense  as  the  body.  This  allows  the  flow  field  to  be  constructed 
utilizing  (1)  stream  tube  integration  of  the  continuity  equation,  (2)  simplification  of  the 
transverse  momentum  equation  which  is  reduced  to  the  stream  normal  momentum  equation,  and 
(3)  oblique  shock  relations. 

For  the  most  part,  those  who  have  used  thin  shock  layer  approximations  have  approached 
the  problem  using  the  indirect  method.  The  technique  described  by  Maslen  (Reference  2-4)  is 
typical.  With  this  method  the  entire  shock  shape  is  assumed  a priori.  Subsequently,  the 
oblique  shock  relations  are  used  to  obtain  the  properties  downstream  of  the  shock.  Then 
each  stream  tube  is  integrated  from  the  shock  to  the  body  utilizing  the  stream  normal  pres- 
sure gradient  equation.  At  some  point  in  this  integration,  a Fictitious  body  which  corres- 
ponds to  the  assumed  shock  shape  is  determined.  In  general,  this  bi.dy  will  not  correspond 
to  the  actual  body  shape.  The  procedure  is  then  to  modify  the  shock  shape  and  repeat  the 
calculations  successively  until  the  computed  body  coincides  with  the  actual  body. 

In  the  present  instance  a somewhat  different  procedure  is  used  under  the  basic  premise 
that  thin  shock  layer  theory  applies.  It  has  been  shown  that  current  techniques  used  to  de- 
termine body  pressure  distributions  perform  adequately  (Reference  2-8).  Rather  than  starting 
with  an  assumed  shock  shape,  the  surface  pressure  ratio,  presented  as  the  ratio  of  local  sta- 
tic pressure  to  stagnation  point  pressure,  is  assumed  known.  Then  with  the  known  pressure 
distribution  at  the  body  surface  and  by  utilizing  the  continuity  equation,  oblique  shock  re- 
lations, and  some  form  of  the  transverse  momentrum  equation,  the  shock  shape  and  shock  layer 
flow  field  can  be  determined  directly. 


2-6 


2-7 


One  additional  point  needs  to  be  noted.  In  view  of  the  fact  that  the  indirect  thin 
shock  layer  method  developed  by  Maslen  does  not  require  that  the  surface  pressure  distribu- 
tion be  be  given,  one  must  ask  what  piece  of  information  could  be  eliminated  if  the  surface 
pressure  is  known?  There  are  several  possibilities  including  relaxation  of  the  requirement 
that  the  shock  jump  conditions  be  satisfied.  This  is,  in  fact,  the  condition  that  is  relaxed 
in  tie  present  method.  That  is,  closure  to  the  problem  as  solved  here  Js  achieved  by  a con- 
tinuity balance  between  the  free  stream  flow  passing  through  the  shock  and  the  flow  in  the 
shock  layer.  The  local  shock  slopes  obtained  in  this  manner  do  not  guarantee  compatibility 
with  the  jump  relations.  This  formulation  employs  a specified  surface  pressure  distribution 
and  an  expression  for  the  stream-normal  pressure  gradient.  This  set  of  conditions  is  suffi- 
cient to  achieve  a unique  solution  to  the  problem.  On  the  other  hand,  if  the  shock  jump 
conditions  are  used  as  an  additional  means  of  achieving  closure,  one  can  relax  either  the 
specification  of  the  normal  pressure  gradient  or  the  specified  surface  pressure  condition. 
Relaxation  of  the  pressure  gradient  specification  may  lead  to  a valid  shock  shape  solution 
given  only  the  surface  pressure  distribution  and  using  both  satisfaction  of  continuity  and 
the  shock  jump  relations  to  achieve  closure.  Relaxation  of  the  surface  pressure  speci- 
fication alone  while  retaining  the  pressure  gradient  specification  would  lead  to  a problem 
formulation  given  by  "direct"  thin  shock  layer  methods.  This  question  is  addressed  further 
in  the  conclusions.  Section  2.4. 

The  procedure  described  in  this  report  consists  of  breaking  the  shock  layer  region 
into  two  distinct  portions:  the  stagnation  region  and  the  remaining  body  region.  The  cal- 

culation starts  with  a determination  of  the  flow  field  in  the  stagnation  region. 

2.2.2  The  Stagnation  Region 

It  is  well  known  that  thin  shock  layer  theory  does  not  formally  apply  in  the  vicinity 
of  the  stagnation  point  where  the  flow  does  not  conform  to  the  body  shape.  In  fact,  the 
streamline  curvature  is  opposite  to  that  of  the  body.  To  overcome  these  difficulties  and  to 
determine  the  shock  standoff  distance  at  the  stagnation  point,  a control  volume  analysis  of 
the  stagnation  region  is  performed.  The  essence  of  this  calculation  is  to  enclose  the  stag- 
nation region  within  a control  volume  with  surfaces  consisting  of  the  shock  surface,  body 
surface,  and  closure  surface  generated  normal  to  the  body  surface  (c.f.  Figure  2-2).  The 
objective  of  the  calculation  is  to  position  the  shock  so  that  a mass  flow  balance  exists 
between  flow  passing  through  the  shock  and  flow  passing  through  the  closure  surface,  enabling 
one  to  determine  the  shock  standoff  at  the  stagnation  point.  The  procedural  steps  arc: 
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Figure  2-2.  Schematic  of  stagnation  region  calculation  procedure. 


1.  First,  some  point  on  the  body,  point  A,  in  the  vicinity  of  the  stagnation  point 
is  selected  where  the  flow  has  to  a large  extent  turned  along  the  body  (rather 
than  being  directed  towards  the  body  as  in  the  imediate  neighborhood  of  the  stag- 
nation point).  It  has  been  found  that  a good  criterion  for  selecting  this  point 
is  where  the  body  Mach  number  equals  about  0.3. 

2.  The  shock  shape  in  the  stagnation  region  is  nearly  normal  since  the  stagnation 
region  control  volume  shown  in  Figure  2-2  is  small.  However,  in  order  to  increase 
accuracy,  the  shock  shape  is  assumed  to  be  given  by: 


where  the  shock  coordinates  xg  and  r$  and  normalized  by  the  body  radius  of  curva- 
ture at  the  stagnation  point,  Rn  , and 

Bst 

H = 1 + 2fi  (2-2) 

where  A is  a first  approximation  to  the  nondimensicnal  shock  standoff,  which  is  to 
be  determined.  This  correlation  was  developed  by  Falanga  (Reference  2-9)  to  be 
applied  over  the  spherical  nose  region  for  sphere  cone  shapes  at  high  Mach  number. 
In  the  stagnation  region  the  body  has  a local  radius  of  curvature  which  corres- 
ponds to  some  effective  spherical  nose;  therefore  Equation  (2-1)  is  applicable. 

The  correlation  given  by  Equation  (2-1)  was  developed  using  Maslen's  complete 
thin  shock  layer  analysis  on  a spherically  capped  body. 

3.  Next,  some  fixed  radial  increment,  Ar,  in  the  free  stream  is  chosen  for  stream 
tube  integration  along  the  shock.  As  each  stream  tube  crosses  the  shock,  the 
flow  conditions  just  downstream  of  the  shock  are  determined  by  the  oblique  shock 
relations:* 


Note  that  these  equations  apply  to  perfect-ideal  gases  where  the  isentropic  exponent  is 
identical  to  the  ratio  of  specific  heats,  and  is  independent  of  all  state  variables. 
Although  this  simplification  employed  herein  is  believed  to  yield  acceptably  accurate 
shock  data,  this  should  be  evaluated  in  the  future. 
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4.  Since  the  entropy  and  total  pressure  are  constant  along  a given  stream  tube  down- 
stream of  the  shock,  then  at  the  body  normal  location: 


P = P 

2j  % 


Here,  the  j indicates  conditions  along  the  body  normal  and  the  s indicates  condi- 
tions just  downstream  of  the  shock.  Given  the  static  pressure  at  the  body  normal 
location  and  the  normal  pressure  gradient,*  the  local  Mach  number  is  determined 


j | ypt2  J 


5.  Now,  by  performing  a streamtube  mass  balance  between  the  free  stream  and  the  local 
body  normal,  the  incremental  streamtube  spacing  at  the  body  normal  is  found  from: 


2nr.|Ar  / m p— 


pt  = 2Trr-  An  /m  P. 

2J  \ p‘/2,  % 


An  = Ar 


f)  p 

ri  V % 

7°  p \ p. 

% (,ptl 


See  the  discussion  in  Section  2.2.4  for  the  procedure  for  determining  the  surface  normal 
pressure  gradient. 
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where  , r ,,  are  the  average  radii  to  the  annuli  of  interest  and 


(y  + 1 )Mijsin2fls 

y-T 

Y + ! 
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, - 1) 

and  the  massflow  function  is  gi ven  by 
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where  w = massflow  rale 


A - streamlube  area  • 

6.  The  incremental  spacing.  An,  for  each  streemtube  is  thus  determined.  As  the  com- 
putation proceeds  botn  along  the  shock  and  along  the  body  normal,  there  will  axis 
some  concition  at  which  the  radial  coordinates  of  the  shock  and  of  the  body  normal 
are  just  equal.  When  this  occurs,  the  computation  is  complete  since  this  yields 
sufficient  conditions  for  control  volume  closure.  If  desired,  the  calculated 
value  of  A can  be  re-input  into  Equation  (2-2)  and  the  process  repeated.  However, 
experience  has  shown  that  if  the  initial  guess  of  A Is  off  by  as  much  as  a factor 
of  two,  it  will  not  have  much  of  an  affect  on  the  computed  shock  standoff  and  shock 
shape.  Therefore,  in  an  efforv  to  reduce  computing  time  this  iteration  is  cur- 
rently not  performed. 

2.2.3  Remaining  Body  Region 

The  body  normal  which  forms  the  upper  surface  of  the  stagnation  region  control  volume 
provides  the  starting  line  for  the  remainder  of  the  flow  computations.  Figure  2-3  illustrates 
the  body  region  flow  computational  scheme.  As  shown,  the  sequence  of  the  calculations  pro- 
ceeds in  the  streamwise  direction  from  one  body  normal,  k-1,  t.o  t.he  next  body  normal,  k.  At 
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each  body  normal  location  a streamtuhe  integration  is  performed  from  the  body  to  the  shock. 

The  flow  at  each  normal  is  broadly  classified  into  two  categories:  "old  flow"  and  "new  flow". 
The  term  "old  flow"  refers  to  flow  which  has  already  passed  through  the  shock  and  has  been 
accounted  for  at  the  previous  body  normal.  The  old  flow  computation  can  be  obtained  from 


An 
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(2-12) 

j = constant 


where  k is  the  body  normal  index.  Since  the  mass  flow  function  is  dependent  upon  the  local 
Mach  number,  the  Mach  number  must  be  computed  from  the  local  value  of  static  to  total  pres- 
sure as  given  by  Equation  (2-6),  where 


with  w = wall 

st  * stagnation  point 
1 = free  Stream 

j = index  along  normal 
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and 


.j  + V, 

P.  = P + Vs  An  (2-14) 

j w 3n  / „ ' ' 

J = 1 

where  an  appropriate  expression  for  the  pressure  gradient  must  be  supplied  as  discussed  in 
Section  2.2,4.  In  Equation  (2-14),  the  total  distance  from  the  body  surface  is  taken  as 
being  the  summation  of  all  streamtube  widths  up  to  the  j-1  location  for  the  purpose  of  evalu- 
ating the  local  pressure. 

The  "new  flow"  computation  is  performed  after  the  "old  flow"  computation  has  been  com- 
pleted. The  term  "new  flow"  refers  to  the  fact  that  this  calculation  accounts  for  the  flow 
that  has  just  passed  through  the  oblique  shock.  The  calculation  sequence  is  as  follows: 
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1. 


The  static  pressure  jump  across  the  shock  at  the  k—  normal  is  computed  from  the 
relation 


where  the  first  term  on  the  right  hand  side  is  the  given  wall  pressure  ratio,  the 
second  term  is  evaluated  at  the  stagnation  point  and  the  third  term  is  obtained 
from  Equation  (2-14)  by  using  a suitable  expression  for  the  pressure  gradient 
across  the  shock  layer. 

2.  From  the  static  pressure  ratios  at  k and  k-1  two  temporary  values  of  shock  angle 
are  found 
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m = k,k  - 1 


3.  From  these  temporary  values  of  shock  angle  the  total  pressure  ratio  acrors  the 
oblique  shock  can  be  computed  from 
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and  the  average  total  pressure  ratio  for  the  "new  flow'  is  obtained  from 


(2-18) 

j = NF  = new  flow 


4. 


Along  the  k—  normal  over  the  "new  flow"  interval  the  average  Mach  number  is  ob 
tained  from 


j NF 
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and 


(2-20) 
j = NF 


5.  Tlie  spacing  interval  along  the  normal  for  the  new  flow  is  calculated  from 


6.  Completion  of  the  calculation  is  obtained  on  the  basis  of  continuity  by  matching 
the  flow  crossing  the  shock  to  the  flow  in  the  "new  flow"  increment  in  the  shock 
layer.  That  is,  the  spacing  increment.  An,  which  was  obtained  from  Step  5 is  addeo 
to  the  normal  coordinate  at  the  k—  location.  The  radial  coordinate  of  the  resul- 
tant tip  of  the  normal  is  then  set  equal  to  the  free  stream  radial  coordinate. 

This  equality  is  defired  explicitly  and  therefore  requires  no  iteration.  This 
causes  the  normal  to  close  on  the  shock.  The  expression  for  th<  normal  pressure 
gradient  also  includes  such  a radial  dependence,  r$ , (see  Equation  (2-23));  however, 
for  the  sake  of  computational  speed  the  radial  coordinate  at  the  last  "old  flow" 
streamtube  is  used  instead  of  the  actual  shock  coordinate.  Once  closure  has  occur- 
red on  the  basis  of  the  continuity  balance,  the  shock  slope  in  the  interval  is  ob- 
tained from  a straight  line  fit  between  the  end  points  of  the  k-1  and  k body  nor- 
mals. The  angle  of  this  shock  segment  may  not  coincide  with  Equation  (2-16). 

Forcing  these  angles  to  coincide  would  overspecify  the  problem  for  the  condition 
of  a known  surface  pressure  distribution  (c.f.  the  discussion  In  Section  2,4). 

2.2.4  Normal  Pressure  Gradient  Relation 

It  was  mentioned  in  Section  2.1  that  in  the  original  formulation  of  the  thin  shock 
layer  equations  It  was  assumed  that  the  pressure  variation  across  the  shock  layer  could  be 
expressed  by 
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In  the  present  Instance  this  would  Imply  that  the  pressure  would  always  Increase  from  the 
body  to  the  shock  if  the  flow  curvature  is  based  upon  the  local  body  curvature.  However, 
even  for  the  case  of  a simple  sphere/cone  configuration,  examination  of  exact  results  indi- 
cates that  indeed  the  pressure  does  not  always  increase  from  the  body  to  the  shock.  In  fact, 
in  two  regions  the  opposite  is  true.  In  the  stagnation  point  regioi  the  streamline  curvature 
is  opposite  that  of  the  body  and,  therefore,  the  pressure  decreases  from  the  body  to  the 
shock.  This  is  also  true  along  the  cone  where  the  flow  first  undergoes  compression  by  the 
shock  and  then  a farther  turning  compression  and,  hence,  upward  curvature  through  the  shock 
layer.  This  also  results  in  a static  pressure  decrease  from  the  body  to  the  shock. 

in  a more  recent  report  (Reference  2-5)  Maslen  presented  an  expression  for  the  stream- 
no  niv  I pressure  gradient  which  includes  the  effect  of  a substantial  velocity  component  normal 
to  ihe  body.  If  Maslen' s more  recent  expression  for  the  pressure  gradient  is  evaluated  at 
the  shock,  and  if  the  gradient  is  assumed  constant  over  the  shock  layer,  the  following  result 
is  obtained 
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Equation  (2-23)  represents  the  average  pressure  gradient  from  the  body  to  the  shock  over  the 
shock  layer  thickness,  AN.  In  order  to  avoid  further  iteration  and  to  decrease  computation 
time,  the  shock  layer  thickness,  All,  is  evaluated  at  the  k-1  position.  The  radius  of  curva- 
ture of  the  shock  is  given  by  adding  the  local  body  radius  of  curvature  to  AN.  In  order  to 
check  the  accuracy  of  this  expression,  it  was  compared  with  the  exact  results  for  a 30°  sphere/ 
cone  ar.  - 5.  Figure  2-4  shows  a comparison  of  the  exact  pressure  variation  to  that  com- 
puted by  the  above  relation.  For  the  purpose  to  which  the  pressure  gradient  is  used  in  the 
present  method  these  results  look  very  good.  The  fact  that  the  above  expression  predicts  the 
reversal  of  pressure  trend  is  encouraging. 


2.2.G  Difficulties  Encountered  Near  the  Stagnation  Region 

In  the  course  of  developing  the  method  described  here,  other  pressure  variations 
across  the  shock  layer  were  considered  to  seek  simplification  of  the  problem.  If  one  assumes 
constant  pressure  through  the  shock  layer,  certain  difficulties  sometimes  arise  making  it 
necessary  to  Introduce  a pressure  reduction  parameter  to  continue  the  computation.  Recall 
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that  when  making  the  new  flow  computation  the  temporary  value  of  shock  angle  based  upon  the 
static  pressure  ratio  across  the  shock  Is  given  by  equation 
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From  tnis  equation  it  is  clear  that  there  is  a value  of  P2s/Pl  above  which  the  quantity  under 
the  radical  is  greater  than  unity  if  the  pressure  is  assumed  constant  across  the  shock  layer. 
The  value  P2s/P^  may  frequently  be  above  this  critical  value  in  regions  of  low  velocity.  As 
mentioned,  near  the  stagnation  point  the  flow  curvature  is  likely  to  be  opposite  to  the  body 
curvature;  hence,  the  static  pressure  should  decrease  from  the  body  to  the  shock.  Therefore, 
if  one  assumes  a normal  pressure  variation  different  from  Equation  (2-23)  the  critical  value 
may  be  exceeded  and  the  computation  will  cease.  Thus,  near  the  stagnation  region  the  ability 
to  continue  the  calculation  is  very  sensitive  to  the  surface  pressure  distribution  and  the 
stream-normal  pressure  gradient.  Once  away  from  the  stagnation  region,  this  problem  is  alle- 
viated. In  order  to  circumvent  this  computational  difficulty,  for  uniform  pressure  across 
the  shock  layer,  P2s/P ^ = Pw/Pj  is  multiplied  by  a factor,  PFAC,  where 


PFAC  = 5^. 
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-n—  to  a value  from  which  the  shock  angle,  0 can  be  calculated.  This  expression 
K1 


This  reduces 

is  exact  at  the  stagnation  point  and  serves  or  a good  approximation  near  the  stagnation  point. 


2.2.6  Flow  Angularity  Effects 

In  its  simplest  form,  the  present  method  assumes  streamlines  parallel  to  the  body.  A 
computational  option  exists,  however,  which  accounts  for  flow  angularity  effects.  The  proce- 
dure followed  to  determine  the  variation  of  flow  angle  across  the  shock  layer  is: 

• Calculate  the  flow  deflection  angle,  Ss,  just  behind  the  shock  (see  Figure  2-5) 
with  the  equation 
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a)  Definition  of  flow  angle,  a. 


b)  Effect  of  flow  angularity  on  streamtube  flow  area. 


Figure  2-5.  Representation  of  flow  angularity  within 
shock  layer. 
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where  e$  = shock  angle 

and  i = flow  deflection  angle  at  the  shock. 

• Define  an  angle,  a = 6g  - 6,  (2-27) 

where  6 = local  flow  deflection  relative  to  the  body. 

• Vary  the  effect  of  the  flow  angle,  a,  on  streamtube  flow  area  (from  zero  at  the 
body  to  its  full  value  at  the  shock)  according  to  the  following  relation: 
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where  At1i  + 1 = l°ca^  stream  tube  flow  area  increment  corrected  for  flow 
angularity 

Anj  = local  stream  tube  flow  area  before  flow  angularity  correction 
i = iteration  index 

w"  = fraction  of  total  distance  across  shock  layer  measured  from 
body. 


The  above  expression,  which  makes  the  flow  area  vary  with  the  cosine  of  a,  was  chosen 
on  intuitive  grounds;  however,  it  nearly  represents  the  Tlow  angle  variation  for  a 30°  sphere/ 
cone  at  M = 5 as  shown  in  Figure  2-6  and  one  anticipates  comarably  good  results  for  other 
configurations  and  Mach  numbers.  The  iteration  index,  i,  is  currently  set  at  a maximum  of 
two  iterations  for  the  sake  of  saving  computation  time. 


2.3  RESULTS 

Shock  shapes  calculated  by  the  thin  shock  layer  technique  were  compared  with  the  exact 
shock  solution  for  a sphere/cone  and  with  PANT  Series  B data  for  the  45°  simple  biconic,  con- 
vex biconic  and  triconic  models  (c. f.  Reference  2-8).  Data  from  PANT  Series  J tests  were 
used  for  comparison  for  the  60°  simple  biconic.  Initial  calculations  were  generated  with 
exact  theoretical  and/c-  experiment  1 pressure  distributions  in  order  to  assess  the  validity 
of  the  thin  shock  layer  approach  for  determining  shock  shapes  alone.  Table  2-1  summarizes 
the  different  cases  considered  for  validation  of  the  technique.  Comparisons  include  shock 
calculations  for  which  the  pressure  across  the  shock  layer  is  assumed  constant  and  for  which 
pressure  gradient  is  assumed  constant  and  Is  calculated  as  described  in  Section  2.2,4.  Also 
included  are  calculations  which  allow  for  flow  angle  variation  between  the  shock  and  body. 
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FABLE  2-1 . SUMMARY  OF  SHOCK  SHAPE  COMPARISONS 


Section  2.3.1  briefly  discusses  the  method  used  previous  to  the  one  described  in  this 
report.  Comparisons  of  the  previous  method  to  exact  results  are  given.  These  previous  re- 
sults provide  a basis  upon  which  the  current  method  can  be  compared.  Section  2.3.2  describes 
results  considering  constant  pressure  normal  to  the  surface  and  parallel  flow.  Section  2.3.1 
describes  results  considering  normal  pressure  gradient,  zero  flow  angle  effects  and  Section 
2.3.4  includes  results  considering  pressure  gradient  and  variable  flow  angle  effects.  The 
results  shown  in  these  sections  were  computed  using  the  exact  (theoretical)  pressure  distri- 
bution along  the  body.  Section  2.3.5  describes  results  using  the  PANT  pressure  correlation 
distribution. 

2.3.1  Results  Using  Previous  Method 

Before  discussing  the  results  obtained  with  the  method  developed  in  this  report, 
termed  the  present  method,  it  is  worth  reviewing  those  results  obtained  with  the  previous 
method.  As  discussed  in  Section  2.1  the  "previous"  method  utilized  a simple  correlation  relat- 
ing shock  shape  to  local  body  shape  in  the  subsonic  region  and  a simple  shock  expansion  tech- 
nique in  the  supersonic  region.  Figures  2-7  and  2-8  present  representative  convex/biconic  and 
triconic  results  obtained  with  the  previous  method  compared  to  exact  results  and  the  original 
PANT  shock  shape  correlation.  (As  a reminder,  the  "original"  PANT  shock  shape  correlation 
employed  only  a correlation  between  shock  shape  and  body  shape.)  It  can  be  seen  that  the 
"previous"  method  does  a fairly  good  job  of  predicting  the  shock  shape  for  the  two  cases 
shown.  These  results  are  not  as  good  as  those  obtained  with  the  new  technique  to  be  discus- 
sed in  the  following  section.  In  spite  of  the  apparent  good  results,  it  should  be  recalled 
that  in  order  to  obtain  the  results  shown  in  Figures  2-7  and  2-8,  it  was  necessary  to  deter- 
mine the  stagnation  point  shock  standoff  independently.  Also,  the  previous  method  could  not 
adequately  handle  shapes  that  are  very  blunted  and  exhibit  large  subsonic  regions.  Finally, 
it  should  be  recalled  that  the  previous  method  relies  basically  on  a curve  tracing  scheme 
(c.f.  Reference  2-2).  That  is,  the  shock  shape  is  started  at  the  stagnation  point  and  built 
up  by  successively  adding  line  segments  whose  slopes  are  computed.  As  small  errors  are  intro- 
duced into  such  a method  the  potential  for  building  large  accumulated  errors  is  serious. 

Hence,  in  practice,  when  unusual  body  shapes  may  occur  the  potential  for  large  inaccuracies 
is  great. 

The  method  described  in  this  report  does  not  rely  on  a curve  tracing  scheme.  Instead 
the  technique  is  basically  an  integration  scheme  accounting  for  total  flow  within  the  shock 
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2-7.  Comparison  of  previous  shock  shape  methods 
with  exact  solutions  of  convex  biconic,  M 


of  previous  shock  shape  methods 
solutions  of  triconic,  M = 5. 


layer.  Therefore,  the  present  method  has  the  advantage  of  any  integration  scheme,  namely, 
that  errors  tend  to  be  compensated  out  in  the  integration  process.  This  will  be  discussed 
further  in  the  paragraphs  that  follow. 

2.3.2  Constant  Pressure  - Uniform  flow  Angle  Results 

Figure  2-9a  presents  a comparison  of  the  thin  shock  layer  prediction  with  the  exact 
shock  solution  for  a 30°  sphere/cone  (M^  = 5,  y = 1.4)  under  the  assumptions  of  2ero  normal 
pressure  gradient  and  parallel  flow.  Note  that  the  terminal  shock  location  is  somewhat  dis- 
placed from  the  true  location  and  that  the  shock  bends  sharply  near  the  shoulder.  This  be- 
havior is  probably  explained  by  the  omission  of  the  normal  pressure  gradient.  Figure  2-9b 
presents  a comparison  of  the  thin  shock  layer  prediction  with  PANT  Series  B data  for  a tri- 
conic nosetip  configuration.  Hence,  the  only  serious  deviation  from  the  true  shock  shape 
occurs  in  the  region  of  high  body  curvature  where  the  pressure  gradient  effect  is  impor- 
tant. In  both  cases,  the  constant  normal  pressure,  uniform  flow  angle  assumptions  yield  a 
shock  shape  reasonably  consistent  with  the  data.  The  stagnation  point  stand-off  calculated 
for  the  jphere/cone  is  in  excellent  agreement  while  that  for  the  triconic  is  about  15  percent 
too  low.  It  appears  also  that  in  a large  region  of  high  interest  (that  which  affects  en- 
tropy swallowing  along  the  body),  the  constant  pressure  uniform  flow  angle  assumption  is  a 
good  first  approximation  to  the  problem. 

2.3.3  Pressure  Gradient  - Uniform  Flow  Angle  Results 

Figure  2-10  presents  results  including  the  effect  of  considering  an  internally  computed 
normal  pressure  gradient  across  the  shock  layer  (Equation  (2-23))  but  still  using  the  uniform 
flow  angle  assumption.  Figures  2-10a,  2-10b,  and  2-10c  present  the  results  for  the  sphere/ 
cone,  triconic  and  convex  biconic,  respectively.  Figures  2-1 Od  and  2-10e  present  the  re- 
sults for  the  45°  and  60°  simple  biconics. 

It  can  be  seer,  that  the  inclusion  of  the  normal  pressure  gradient  term  improves  shock 
shapes  calculated  for  the  sphere/cone  and  triconic.  It  is  significant  that  the  improvement 
occurs  mainly  in  regions  of  high  curvature  where  the  absence  of  a normal  pressure  gradient 
causes  kinks  to  develop.  The  stagnation  point  stand-off  for  the  tricoric  has  actually  been 
further  decreased  a small  amount  by  the  inclusion  of  a normal  pressure  gradient  term. 

The  convex  biconic  and  the  45"  simple  biconic  shock  calculations  show  excellent  agree- 
ment with  the  exact  calculation,  The  60°  biconic  calculation  shows  somewhat  erratic  behavior. 
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Figure  2-9.  Results  of  thin  shock  layer  calculation  assuming 
constant  pressure  - zero  flow  angularity  (M  = 5, 
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Figure  2-10.  Concluded 


Since  the  flow  along  the  60°  biconic  is  subsonic  along  the  entire  60°  face,  the  whole  region 
could  be  considered  a stagnation-like  region.  Hence,  it  is  expected  that  nonuniformities  in 
flow  direction  are  present  all  along  the  60°  face  since  the  shock  remains  strong  to  the  shoul- 
der. Thus,  it  is  suggested  that  not  considering  the  flow  angle  variation  accounts,  in  part, 
for  the  low  stagnation  point  stand-off  prediction.  In  turn,  the  underprediction  of  the  stagna- 
tion point  shock  stand-off  may  account  for  the  oscillatory  behavior  away  from  the  stagnation 
point.  This  oscillation  is  hypothesized  to  result  from  the  tendency  of  the  flow  to  account 
for  the  correct  total  mass  flow  between  the  shock  and  the  body.  The  possibility  also  exists 
that  this  oscillatory  behavior  may  also  be  the  result  of  some  numerical  problem  resulting  from 
spacing  considerations  of  the  local  body  normals  or  from  implementation  of  the  closure  model, 
as  discussed  in  Sections  2.2.1  and  2.2.3.  Preliminary  studies  of  nodal  spacing  were  performed. 
The  results  of  this  study  did  not  conclusively  indicate  that  nodal  spacing  is  tne  cause  of  the 
oscillatory  behavior,  further  studies  are  required  to  resolve  this  uncertainty. 

2.3.4  Pressure  Gradient  - How  Angle  Variation  Results 

figure  2-11  presents  results  which  include  the  effect  of  a finite  pressure  gradient  as 
discussed  in  Section  2.2.4  and  flow  angularity  effects  across  the  shock  layer  as  discussed  in 
Section  2.2.5.  figure  2- 1 1 a and  2-llb  present  these  results  for  the  sphere/cone  and  triconic. 
Figure  2-1 1 c and  2-lld  present  the  results  for  the  45°  and  60°  simple  biconics.  The  accuracy 
of  the  results  obtained  by  including  the  flow  angle  variation  shows  no  substantial  improvement 
and  perhaps  somewhat  poorer  results  than  those  discussed  in  Section  2.3.3.  Note  that  in  Fig- 
ure 2-lld  there  again  is  considerable  waviness  in  the  shock  for  the  60°  simple  biconic.  Again, 
the  assumption  is  that  due  to  the  fact  that  the  stagnation  point  shock  stand-off  is  considerably 
underpredicted,  the  shock  shape  tends  to  compensate  towards  the  true  slope  in  a somewhat 
oscillatory  manner. 

Comparing  the  results  in  detail  shown  in  Figure  2-10  with  those  shown  in  Figure  2-11, 
it  can  be  seen  that  the  inclusion  of  the  flow  angle  variation  has  little  effect  on  the  overall 
prediction.  In  the  region  of  the  stagnation  point  the  flow  angularity  apparently  tends  to 
counteract  the  negative  pressure  gradient  with  regard  to  shock  stand-off.  This  explains  the 
good  agreement  between  the  present  method  and  the  exact  shock  solution  attained  in  Figure  2-9 
by  assuming  constant  pressure  and  zero  flow  angularity  across  the  shock  layer. 
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2.3.5  Results  Using  PANT  Correlation  Pressure  Distributions 


Because  the  major  objective  of  this  task  is  to  determine  the  applicability  of  the  thin 
shock  layer  method  to  the  PANT  shape  change  code,  several  shock  calculations  were  made  using 
PANT  correlation  generated  pressure  distributions.  Hence,  these  results  indicate  the  overall 
agreement  obtained  by  use  of  the  thin  shock  layer  method  in  conjunction  with  pressure  distri- 
butions determined  by  correlative  methods.  Figure  2-12  presents  results  using  the  PANT  corre- 
lation pressure  distribution  in  place  of  the  exact  pressure  distributions  (pressure  gradient 
and  flow  angularity  effects  were  included  for  these  calculations).  Figure  2-13  presents  com- 
parisons of  the  exact  and  PANT  correlation  pressure  distributions  for  each  body.  Figure  2-12a, 

2- 1 2b , and  2-12c  present  the  results  for  the  triconic,  45°  simple  biconic  and  60°  simple  biconic, 
respectively.  The  comparisons  show  that  for  the  most  part  the  PANT  pressure  correlation  is 
adequate  for  use  in  the  thin  shock  layer  method.  However,  it  should  be  noted  that  in  Figure 
2- 1 2a  and  2-12.:  some  kinks  in  the  shock  shape  developed.  These  were  traced  back  to  body  point 
spacing  in  relation  to  the  smoothness  of  the  pressure  correlation  distributions.  This  effect 
deserves  further  consideration  before  incorporating  into  the  PANT  code.  The  60°  biconic  shock 
calculation  still  does  not  compare  with  the  data  as  favorably  as  do  the  calculations  for  the 
other  geometries,  but  in  a gross  sense  it  remains  an  improvement  over  tn<  present  PANT  shock 
determination  techniques  as  discussed  in  Section  2.3.1. 

Since  shock  angle  and  not  shock  shape  are  of  primary  importance  in  determining  boundary 
layer  edge  conditions,  the  shock  angle  variation  was  investigated.  Figure  2-14  presents  plots 
of  shock  angle  versus  shock  radius  for  the  45°  simple  biconic  and  triconic  using  the  PANI  sur- 
face pressure  correlations.  Each  plot  compares  the  exact  data  with  prediction  using  the  cur- 
rent method.  These  plots  give  a somewhat  better  indication  of  the  overall  effect  upon  the 
entropy  swallowing  calculation. 

In  each  case,  the  comparison  indicates  a slight  dip  in  the  predicted  shock  angles 
located  in  the  vicinity  of  the  stagnation  region.  This  dip  corresponds  to  the  negative  shoct 
curvature  found  near  the  nose  in  Figures  2-I2a  and  2-1 2b  and  seems  to  correspond  directly  to  a 
slight  underprediction  of  the  stagnation  point  shock  stand-off.  It  is  interesting  to  note, 
however,  that  the  method  attempts  to  compensate  for  this  Initial  mismatch. 

2.4  CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  obtained  with  this  method  arr-  generally  good.  In  a gross  sense  these  results 
are  better  than  those  obtained  with  previous  metnods  and  the  method  appears  to  offer  a powerful 
approach  to  the  prediction  of  flow  in  the  shock  layer.  The  method  is  also  numerically  efficient 
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Figure  2-14.  Variation  of  shock  angle  with  radius  using 
PANT  pressure  correlations  (M  - 5,  , = 1.4). 
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and  requires  significantly  less  computing  time  than  highly  sophisticated  procedures  without 
sacrificing  overall  accuracy  or  generality.  A typical  computation  involving  AO  integration 
normals  takes  about  2 seconds  of  computer  running  time  on  the  Univac  1108.  Perhaps  the  most 
desirable  aspect  of  the  method  described  in  this  report  is  that  covergence  of  a solution  is 
assured.  That  is,  a shock  shape  prediction  method  which  is  used  in  a shape  change  code  must 
be  utilized  many  tens  of  times  during  a trajectory.  A method  for  which  convergence  is  sensi- 
tive to  body  shape  (for  example)  would  not  serve  the  desired  purpose. 

As  noted,  when  the  shock  stand-off  is  underestimated  for  blunt  shapes  a considerable 
amount  of  waviness  in  the  shock  may  occur.  Such  shock  waviness  may  adversely  affect  the  bound- 
ary layer  calculations.  It  is  felt  that  the  basic  cause  of  the  shock  waviness  can  be  resolved. 

A part  of  the  problem  may  be  related  to  the  explicit  nature  of  some  of  the  calculational  pro- 
cedures relating  to  normal  pressure  gradient  and  flow  angularity.  It  should  be  recalled  that 
for  the  sake  of  computational  speed  certain  assumptions  were  made  regarding  the  necessity  to 
update  variables  relating  to  normal  pressure  gradient  and  flow  angularity.  These  assumptions 
would  lead  to  some  error  accumulation  in  the  system  that  will  propagate  through  the  flow.  These 
effects  may  be  curtailed  by  reducing  step  size  or  by  upgrading  the  iterative  procedures. 

It  was  also  mentioned  that  at  a sharp  shoulder  point  where  the  combined  effects  of  high 
Mach  number  and  large  curvature  are  both  present  causes  problems  attributable  to  the  inadequacy 
of  the  normal  pressure  gradient  expression.  The  problem  relating  to  definition  of  the 
normal  pressure  gradient  in  regions  of  large  curvature  and  Mach  number  coincides  with  the 
fact  that  stream  curvature  effects  are  actually  propagated  along  Mach  lines  from  the  body  to 
the  shock.  Consequently,  in  these  regions  it  may  be  more  beneficial  to  perform  the  streamtube 
integrations  along  the  Mach  lines  rather  than  along  body  normals.  This  would  more  properly 
describe  the  physical  behavior  of  the  flow  within  the  shock  layer. 

It  was  mentioned  earlier  in  Sections  2.2.1  and  2.2.3  that  if  the  surface  pressure  dis- 
tribution is  given  that  one  of  the  conditions  for  closure  to  the  problem  can  be  relaxed.  In 
the  present  method  the  requirement  of  totally  satisfying  the  oblique  shock  jump  condition  is 
relaxed.  If  instead,  the  surface  pressure  distribution  obtained  from  PANT  correlations  was 
only  used  as  a first  approximation  to  the  actual  surface  pressure  it  would  be  possible  to 
correct  the  surface  pressure  by  demanding  closure  based  upon  both  continuity  and  satisfaction 
of  the  oblique  shock  relations.  The  surface  pressure  distribution  could  be  corrected  in  a step 
by  step  procedure  from  the  stagnation  point  to  the  sonic  line.  Additionally,  a closure  condi- 
tion based  upon  the  elliptic  nature  of  the  subsonic  region  would  have  to  be  imposed.  One 
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possible  condition  for  this  closure  would  be  the  requirement  of  satisfying  global  momentum 
over  this  region. 

A second  approach  to  the  problem  was  also  mentioned  in  Section  2.2.1.  Namely,  if  the 
surface  pressure  distribution  is  known  closure  may  be  based  upon  the  satisfaction  of  continuity 
and  oblique  shock  relations  if  the  specification  of  the  normal  pressure  gradient  is  eliminated. 
The  normal  pressure  would  be  assumed  to  vary  linearly  but  with  no  predefined  slope.  Essen- 
tially, satisfaction  of  the  closure  conditions  will  determine  the  value  of  the  normal  pressure 
gradient.  Solution  convergence  of  such  a method  is  not  known  at  present. 

As  a result  of  this  study  several  recommendations  seem  in  order: 

1.  The  problems  of  shock  waviness  and  normal  pressure  gradients  near  shoulder-like 
regions  require  further  investigation.  Suggestions  to  correct  these  problems 
were  made  in  the  previous  paragraphs.  It  is  recommended  that  these  procedures 
be  adopted. 

2.  An  investigation  of  the  effect  of  shock  shape  and  body  pressure  distribution  wavines 
and  accuracy  upon  resultant  heat  transfer  and  shape  change  calculations  should  be 
performed.  The  results  of  such  a sensitivity  study  will  show  whether  small  amounts 
of  shock  waviness  can  be  accepted.  It  would  also  indicate  the  combined  accuracy 
requirements  for  shock  shape  and  boundary  layer  edge  conditions.  Such  a sensitivity 
study  should  also  investigate  the  effects  of  nodal  spacing  and  the  implications  upon 
the  surface  pressure  correlation  smoothness. 

3.  It  was  mentioned  in  this  section  that  another  approach  to  the  problem  would  be  to 
eliminate  specification  of  the  normal  pressure  gradient  and  substitute  instead  a 
closure  condition  based  upon  satisfying  the  oblique  shock  relations.  It  Is  recom- 
mended that  this  approach  be  attempted  in  that  the  oblique  shock  relations  are  known 
more  exactly  than  the  normal  pressure  gradient.  This  technique  could  then  be  used 
as  a standard  for  obtaining  shock  shapes  based  upon  exact  surface  pressure  distri- 
butions. 

4.  It  is  recommended  that  the  other  technique  mentioned  in  this  section,  namely,  the 
method  that  solves  for  the  surface  pressure  distribution  be  developed.  The  results 
obtained  from  the  work  under  recommendation  (3)  will  be  used  as  a standard  of  com- 
parison in  that  the  shock  shapes  are  obtained  from  exact  surface  conditions.  This 
would  be  of  great  aid  in  the  development  of  the  required  normal  pressure  gradient 
expression  required  to  perform  recommendation  (4). 


5.  Investigate  the  techniques  required  to  account  for  a step  change  In  y across  the 


shock.  Modify  the  calculations  to  account  for  these  real  gas  effects. 
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SECTION  3 

BOUNDARY  LAYER  TRANSITION  ON  SMOOTH  NOSET  IPS 
by 
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LIST  OF  SYMBOLS 


dimensionless  blowing  parameter.  Equation  (3-7) 


heat  transfer  coefficient 


specific  turbulent  energy 


frequency 


turbulence  intensity,  I - Vu'Vu 


roughness  height 


Mach  number 


momentum  thickness  Reynolds  number 


spheriuai  nose  radius 


wetted  distance  from  stagnation  point 


absolute  temperature 


turbulent  velocity  fluctuation 


streamwise  velocity 


velocity  normal  to  surface 


missile  velocity 


ballistic  coefficient 


entry  angle 


combined  disturbance  parameter.  Equation  (3-10) 


kinematic  viscosity 


specific  turbulent  dissipation  rate 


wall  cooling  ratio  parameter,  Equation  (3-6) 
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LIST  OF  SYMBOLS  (CONCLUDED) 


P 


den'  i ty 


e momentum  thickness 


Subscripts 

e evaluated  at  boundary  layer  edge 

T transition  point 

w wall 
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3.1 


INTRODUCTION 


Boundary  layer  transition  on  passive  missile  nosetips  has  important  implications  to 
mission  accomplishment.  The  resulting  transitional/turbulent  boundary  layer  causes  increased 
surface  recession  and  surface  temperatures  which  require  adequate  ablative  material  to  pro- 
tect the  payload.  The  greater  the  amount  of  ablative  material  requried  and  the  higher  the 
surface  temperature,  the  more  likely  the  tips  will  suffer  thermostructural  failure. 

Nosetip  boundary  layer  transition  initiated  by  roughness  induced  disturbances  is  reason- 
ably well  understood  (see  Reference  3-1).  However,  the  quantification  of  distrubances  arising 
from  other  sources  in  the  reentry  environment  and  their  influence  on  nosetip  transition  is  not 
well  understood.  Higher  strength  graphites  now  under  development  have  a finer  microstructure 
than,  for  example,  ATJ-S  and  therefore  evolve  smoother  surfaces  during  laminar  ablation.  Transi- 
tion on  nosetips  fabricated  from  these  smoother  graphites  may  very  well  not  be  roughenss  domin- 
ated so  that  the  performance  of  such  nosetips  is  somewhat  uncertain. 

An  example  prediction  of  roughness  dominated  nosetip  transition  is  presented  in  Fig- 
ure 3-1  which  shows  the  transition  onset  altitude  as  a function  of  graphite  surface  roughness 
for  a typical  ballistic  reentry  nosetip  and  trajectory.  Nosetip  transition  onset  refers  to  the 
first  occurrence  of  tu-bulent  boundary  layer  flow  ahead  of  the  sonic  point,  an  event  which  is 
crucially  important  to  nosetip  recession  and  thermostructural  performance.  Current  graphites, 
ATJ-S  for  example,  have  peak-to-val ley  roughness  heights  of  about  0.5  mil  which  result  in 
typical  transition  altitudes  in  the  order  of  50,000  feet  or  so.  For  smoother  surfaces,  transi- 
tion occurs  at  lower  altitudes.  It  is  more  difficult  to  predict  transition  altitude  for  smoother 
materials  because  of  the  probable  influence  of  other  disturbances.  The  present  study  was  under- 
taken in  order  to  improve  the  understanding  of  this  "smooth  wall"  nosetip  transition  problem, 

3.1.1  Objectives 

In  order  to  provide  a technique  for  improved  prediction  of  transition  on  finer  micro- 
structure  graphite  nosetips,  the  following  objectives  were  established. 

• To  generalize  the  PANT  rough  wall  transition  criteria  to  include  free  stream  dis- 
turbances 

• To  begin  characterization  of  the  free  stream  disturbance  to  which  nosetips  are  sub- 
jected in  the  atmosphere. 
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fil'jre  3-1.  Roughness  dominated  nosetip  boundary  layer  transition 
onset. 
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The  first  objective  is  stated  in  terms  of  the  well  established  rough  wall  criteria  because  of 
the  convenience  and  desirability  of  using  a single,  continuous  expression  over  the  entire  range 
of  roughness  influence. 

3.1.2  Approach 

The  current  study  consisted  of  a review  of  transition  data  from  relatively  smooth  sur- 
faces, the  extension  of  the  rough  wall  disturbance  parameter  to  include  free  stream  turbulence 
and  a brief  assessment  of  the  reentry  free  stream  disturbance  environment. 

The  review  of  pertinent  transition  data  is  presented  in  Section  3.2.  Transition  data 
taken  in  supersonic  wind  tunnels  is  excluded  from  consideration  because  of  the  unknown  influ- 
ence of  noise  radiated  from  the  tunnel  wall  boundary  layers.  For  smooth  blunt  bodies,  there 
remains  for  consideration  a few  flight  transition  teats,  one  transonic  wind  tunnel  experiment 
and  one  shock  tube  test.  Data  obtained  in  subsonic  wind  tunnels  on  smooth  flat  plates,  for 
which  stream  disturbances  are  relatively  well  documented,  are  also  reviewed. 

Section  3.3  addresses  the  modeling  of  transition  on  blunt  bodies  subject  to  disturbances 
in  the  air  stream  and/or  disturbances  arising  from  surface  roughness.  The  PANT  rough  wall  cri- 
teria are  discussed.  The  results  of  a theoretical  transition  prediction  technique  for  flat 
plates  and  blunt  bodies  with  smooth  surfaces  are  renewed.  An  extension  of  the  rough  wall  dis- 
turbance parameter  to  include  stream  disturbances  ia  suggested.  This  leads  to  a revised  cri- 
terion which  takes  into  account  the  influence  of  stream  turbulence  as  well  as  wall  roughness. 

A brief  assessment  of  the  disturbances  present  in  the  atmosphere  is  presented  in  Section 
3.4.  Consideration  is  given  to  vehicle  vibration,  clear  air  turbulence  and  cirrus  ice  clouds. 

A summary  of  the  results  of  the  present  study  is  given  in  Section  3.5  together  with  the 
conclusions  drawn  from  these  results. 

3.2  SMOOTH  WALL  TRANSITION  DATA 

Boundary  layer  transition  is  the  process  during  which  some  disturbance  is  amplified  into 
turbulence.  Since  the  response  of  the  boundary  layer  may  be  different  to  different  magnitudes 
and  types  of  disturbances,  the  understanding  of  transition  hinges  on  the  understanding  of  the 
existing  disturbances.  It  follows  that  the  usefullness  of  any  transition  data  is  limited  by 
the  extent  to  which  the  disturbance  environment  is  known.  For  smooth  walls  (where  transition 
is  not  roughness  dominated)  the  measurement  of  stream  disturbances  is  important.  It  has  not, 
however , always  been  regarded  as  an  essential  pare  of  transition  experiments.  Therefore,  not 
all  smooth  wall  transition  data  are  considered  to  be  useful  herein. 
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In  supersonic  wind  tunnels,  transition  on  smooth  surfaces  is  usually  dominated  by  sound 
radiated  from  turbulent  tunnel  wall  boundary  layers.  Since  noise  (acoustic  flow  field  fluctua- 
tions) is  not  believed  to  be  an  important  source  of  disturbance  in  the  reentry  environment, 
noise  dominated  transition  data  are  of  little  interest.  For  that  reason,  data  obtained  in 
supersonic  wind  tunnels  have  been  excluded  from  consideration  herein. 

Included  for  consideration  in  the  present  study  are  data  from  several  flight  tests 
(References  3-2,  3-3,  3-4,  3-5,  and  3-6)  and  the  subsonic  shrounded  model  data  in  Dunlap  and 
Kuethe  (Reference  3-7).  The  blunt  body  shock  tube  data  of  Stetson  (Reference  3-8)  are  also 
included  although  stream  turbulence  levels  were  not  documented.  Subsonic  flat  plate  data, 
for  which  turbulence  levels  are  reasonably  veil  documented,  are  also  of  interest  (References 
3-9  and  3-10). 

In  the  analysis  of  the  blunt  body  data,  peak-to-val ley  values  of  the  surface  roughness 
height  were  taken  to  be  4 times  the  reported  RMS  values.  The  justification  for  this  procedure 
is  shown  in  Figure  3-?  where  peak-to-valley  surface  roughnesses  for  the  PANT  rough  wall  calor- 
imeters are  compared  to  RMS  values  measured  with  a profilometer. 

3.2.1  NACA  Calorimeter  Flight  Tests 

These  tests  are  especially  significant  in  that  they  demonstrated  that  very  high  values 
of  the  momentum  thickness  Reynolds  number  at  transition  (Refj  T ~ 1000)  can  be  obtained  on  care- 
fully polished  metal  blunt  bodies  in  a low  disturbance  environment.  These  transition  Reynolds 
numbers  are  about  the  same  as  are  observed  for  smooth  flat  plates  in  subsonic  wind  tunnels 
with  low  stream  turbulence  levels  (these  flat  plate  data  are  discussed  in  Section  3.2.6).  The 
often  quoted  nominal  value  of  Ree  ^ - 300  for  blunt  bodies  obviously  does  not  apply  to  this 
situation. 

Three  flight  tests  performed  by  the  Langley  Aeronautical  Laboratory  of  NACA  are  documented 
in  References  3-2,  3-3,  and  3-4.  The  nosetips  tested  were  inconel  spherical  shell  calorimeters 
instrumented  with  back  wall  thermocouples.  In  two  cases  transition  data  were  obtained  during 
both  the  powered  and  coasting  portions  of  the  trajectory. 

Reference  3-2  reported  a surface  roughness  of  0.025  mil  RMS  as  measured  with  a profilom- 
eter. Surface  roughenss  of  0.005  mil  RMS  and  0.003  mil  RMS  measured  with  an  interference  micro- 
scope were  reported  in  References  3-3  and  3-4,  respectively . It  is  believed  Lliat  profilometer 
measurement  of  relatively  smooth  surfaces  are  inaccurate  because  (1)  the  stylus  is  about  the 
same  size  or-  larger  than  the  surface  irregularities  and  (2)  the  stylus  itself  may  smooth  out 
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irregularities  in  its  path.  Therefore,  the  surface  roughness  stated  in  Reference  3-2  should  be 
regarded  as  a lower  bound  on  the  actual  roughness. 

Representative  transition  data  from  these  flights  are  listed  in  Table  3-1.  The  data 
used  here  are  limited  to  transition  locations  on  the  spherical  nose.  EROS  boundary  layer  solu- 
tions were  used  to  generate  the  momentum  thicknesses  and  Reynolds  numbers.  The  calculated 
values  of  Reg  were  found  to  be  in  reasonable  agreement  with  those  reported  in  References  3-2, 
3-3,  and  3-4.  As  is  explained  in  Section  3.3  (or  Reference  3-1),  the  parameter  (Te/Tw)(k/fi) 
is  a measure  of  the  disturbance  introduced  by  the  surface  roughness. 

For  the  flight  reported  in  Reference  3-2,  the  transition  Reynolds  numbers  are  fairly 
low.  This  is  believed  to  be  primarily  a roughness  influence.  The  Re0  T of  435  obtained  with 
the  engine  off  would  be  predicted,  using  the  rough  wall  criterion,  for  a surface  roughness  of 
0.15  mil.  This  apparent  error  of  30  to  50  percent  is  believed  to  be  well  within  the  combined 
uncertainty  of  the  surface  roughness  measurement,  the  factor  of  4 conversion  to  peak-to-val ley 
value  and  the  transition  criterion  itself.  It  would  require  a surface  roughness  of  0.20  mil 
to  predict  the  Ren  T = 236  observed  with  the  engine  on.  Apparently,  vehicle  vibrations  from 

a i I 

the  engine  enhanced  the  influence  of  roughness  by  a factor  of  4/3. 

The  surface  roughness  reported  in  References  3-3  and  3-4,  if  accurate,  would  be  expected 
to  have  little  influence  on  transition  (predicted  roughness  induced  Re0  T > 3000).  Vehicle 
vibration  (rocket  engine  on)  did  not  cause  particularly  low  Re(j  T in  the  absence  of  significant 
roughness.  No  explanation  is  apparent  for  the  lower  ReQ  j for  the  smoother  surface  reported  in 
Reference  3-4  as  compared  to  Reference  3-3.  The  level  of  the  transition  Reynolds  numbers  for 
these  two  flights  give  some  indication  of  the  magnitude  of  atmospheric  disturbance  below 
10,000  feet. 

For  all  three  of  these  flights,  momentum  thickness  Reynolds  numbers  at  transition  were 
much  larger  when  transition  occurred  on  the  aft  surfaces  rather  than  on  the  spherical  portion 
of  the  model.  A maximum  Re0  j of  about  3500  was  reported  in  Reference  3-3.  These  very  large 
values  of  transition  Reynolds  numbers  may  indicate  that  the  turbulence  intensity  at  the  bound- 
ary layer  edge  along  the  aft  surfaces  (but  not  necessarily  on  the  nose)  was  much  lower  than 
that  ever  achieved  in  wind  tunnel  facilities  (see  Section  3.2.6). 

3.2.2  Lockheed  X-17  Flight  Tests 

Transition  data  from  these  powered  reentry  flight  tests  were  obtained  from  Reference  3-5. 
lwo  flights,  R-2  and  R-9  are  described.  The  two  nosetips  were  copper  hemisphere-cylinder  shells 
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TABLE  3-1.  LANGUV  FLIGHT 


Ref. 

kf 

(mil ) 

rn 

(in. ) 

Altitude 

i 

(103  ft) 

V 

(103  Ik) 

sT/RN 

2 

0.10 

4.0 

3.4 

2.8 

0.26 

2 

0.10 

4.0 

7.2 

3.1 

0.58 

3 

0.020 

4.0 

2.5 

2.6 

1.40 

3 

0.020 

4.0 

6.9 

3.1 

1.57 

4 

0.012 

i 

6.5 

10.0 

3.4 

0.78 

Peak-to-valley  (PTV). 
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DATA 


ft!) 

Re0,T 

Rocket 

Engine 

0.43 

235 

on 

0.24 

435 

off 

0.02 

1,200 

on 

0.01 

1,400 

off 

0.02 

800 

on 

with  nose  radii  of  4.5  inches.  Thermocouples  were  installed  in  copper  plugs  shrunk  fit  into 
holes  in  the  shell  prior  to  surface  finishing.  The  surface  finish  on  R-2  was  polished  copper 
with  k = 0.002  mil  RMS.  The  R-9  nosetip  was  nickel  plated  and  polished  to  obtain  a k = 0.001 
mil  RMS  surface  roughness.  A patch  of  k = 0.030  mil  RMS  surface  roughness  was  installed  over 
one  streamwise  ray  of  thermocouples  on  R-9.  The  method  of  determining  surface  roughness  was  not 
reported . 

A limited  but  representative  sample  of  transition  data  from  these  two  flights  is  listed 
in  Table  3-2.  All  data  are  from  the  coasting  portion  of  the  flight.  Again,  the  EROS  code  was 
used  to  generate  momentum  thicknesses  and  Reynolds  numbers.  No  values  of  Reft  T were  reported 
in  Reference  3-5  for  comparison. 

Transition  Reynolds  numbers  for  flight  R-2  are  quite  low  for  the  surface  finish  stated. 
For  flight  R-9,  transition  on  the  roughened  ray  is  in  good  agreement  with  the  PANT  rough  wall 
criterion  which  predicts  transition  20  percent  closer  to  the  stagnation  point.  On  the  highly 
polished  section,  transition  is  delayed  until  Reg  = 800  which  is  in  reasonable  agreement  with 
the  results  discussed  in  Section  3.2.1. 

The  explanation  for  the  factor  of  2 between  Ree  j on  smooth  surfaces  for  the  two  flights 
may  lie  in  the  method  of  thermocouple  installation.  The  thermocouple  plugs,  when  unplated, 
as  in  flight  R-2,  may  not  have  been  polished  evenly  with  the  rest  of  the  surface.  Since  the 
thermocouples  were  installed  on  streamwise  rays,  unevenly  polished  plug  installations  could 
very  well  affect  transition  at  downstream  thermocouples.  Different  atmospheric  condit'ons  for 
the  two  flights  also  might  provide  an  explanation  for  tP^>  discrepancy  in  the  two  values  for 
the  smooth  wall  ReQ  p 

3.2.3  G.E.  Mark  2 Reentry  Flight  Tests 

Two  flight  tests  of  the  G.E.  Mark  2 reentry  body  are  reported  in  Reference  3-6.  Since 
information  about  trajectories,  body  size,  surface  roughness  and  thermocouple  installation 
were  not  presented  there,  these  flight  tests  will  only  be  discussed  briefly  here. 

These  nosetips  were  solid  copper  heat  sinks  with  50°  cone  half  angle  and  spherical  noses. 
Momentum  thickness  Reynolds  numbers  at  transition  were  reported  to  be  in  the  range  of  100  to 
300  for  both  flights.  Again,  the  explanation  of  such  low  values  of  Ren  T probably  lies  with 

U t I 

unknown  surface  irregularities. 
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TABLE  3-2.  LOCKHEED  X-17  FLIGHT  TEST  DATA 


Flight 

k+ 

(mil ) 

rn 

(in. ) 

Altitude 
(103  ft) 

V 

(">•  Sc) 

Vrn 

R-2 

0.008 

m 

31.0 

10.2 

0.45 

R-9 

0.004 

Hi 

31.0 

9.0 

1.50 

R-9 

0.120 

4.5 

31.0 

9.0 

0.45 

3.2.4  Dunlap  and  Kuetne  Shrouded  Model  Data 

This  experiment,  documented  in  Reference  3-7,  was  conducted  using  a 4.5-inch  radius 
spherical  model  in  a transonic  shroud  thereby  obtaining  hypersonic  pressure  distributions 
without  the  penalty  of  the  high  noise  levels  present  in  hypersonic  tunnels.  The  models  were 
machined  from  brass  and  polished  to  a finish  of  k = 0.001  mil  RMS  as  measured  with  a profilom- 
eter.  Turbulence  intensities  were  reported  to  bo  less  than  0.001,  referred  to  local  air  speed 
in  the  subsonic  region.  The  model  was  internally  cooled  in  order  to  achieve  wall  temperature 
ratios  in  the  range  0.4  < Tw/Tfi  < 1.0.  Local  transition  was  detected  using  a Stanton  tube 
mounted  in  the  sonic  region. 

Transition  was  reported  to  occur  in  the  range,  550  < Re0  y < 700,  independent  of  wall 
cooling  ratio.  If  surface  roughness  were  exerting  a significant  influence,  the  Reynolds  num- 
ber at  transition  would  be  expected  to  decrease  with  decreasing  wall  cooling  ratio.  It  is 
therefore  probable  that  transition  was  initiated  by  disturbances  in  the  air  stream. 

A second  conclusion  can  be  drawn  from  the  independence  (in  this  experiment)  of  Re(j  y 
from  the  wall  cooling  ratio.  Linear  stability  analyses  for  smooth  walls  indicate  that  a lower 
wall  cooling  ratio  should  have  a stabilizing  effect  on  the  transition  Reynolds  number  through  its 
influence  on  boundary  layer  profile  shapes.  This  effect  was  not  observed  in  these  tests.  The 
profile  influence  is  also  predicted  for  wall  mass  transfer.  The  observed  independence  of  Re^  y 
from  the  wall  cooling  racio  therefore  probably  also  implies  independence  from  wall  mass  trans- 
fer, at  least  fo"  moderate  blowing  rates. 

3.2.5  Stetson  Shock  Tube  Data 

This  shock  tube  transition  experiment  is  documented  in  Reference  3-8.  Small  quartz 
hemispheres,  1/4  and  1/2  inch  in  radius,  instrumented  with  thin  film  heat  transfer  gages  were 
tested  at  wall  cooling  ratios  less  than  0.1.  The  surface  roughness  height  W3S  reported  to  be 
less  than  0.001  mil  RMS.  However,  the  existence  of  some  pits  and  scratches  about  0.01  mil 
deep  was  noted.  No  measurements  of  the  shock  tube  turbulence  level  were  reported. 

Transition  occurred  on  the  hemispheres  at  Re0  » 225  for  Tw/Te  * 0.1  increasing  to 
Refl  * 325  for  Tw/Tg  * 0.03.  This  trend  is  opposite  to  that  expected  for  roughness  dominated 
transition.  Shock  tube  conditions  (rather  than  model  temperature)  were  varied  in  order  to 
vary  the  wall  cooling  ratio.  It  seems  likely  that  an  associated  variation  in  turbulence  level 
caused  the  change  in  Re0  y. 
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3.2.6  Low  Speed  Flat  Plate  Data 


Although  these  data  were  not  taken  on  models  of  the  desired  geometry,  they  are  the 
data  for  which  the  best  documentation  of  disturbances  is  available,  In  Section  3.3  it  will 
be  shown  how  these  data  can  be  interpreted  for  application  to  smooth  wall  transition  on  blunt 
body  configurations. 

The  momentum  thickness  Reynolds  number  at  transition  is  shown  as  a furction  of  free 
stream  turbulence  intensity  in  Figure  3-3.  The  earlier  data,  which  represent  the  work  of 
several  investigators  (Schubauer  and  Skramstad,  Hall  and  Hislop,  Dryden),  are  taken  from 
Reference  3-9.  More  recent  data,  from  the  work  of  Spangler  and  Wells,  Reference  3-10,  are 
also  shown. 

The  work  of  Spangler  and  Wells  is  most  complete  in  that  distributions  of  disturbance 
energy  with  frequency  (energy  spectra)  were  measured  and  reported.  Schubauer  and  Skramstad 
also  reported  some  spectral  information.  According  to  linear  stability  theory  (see  Refer- 
ence 3-9,  for  example)  disturbances  in  a certain  frequency  range  are  amplified  while  all 
others  are  damped.  The  energy  within  this  range  is  most  likely  a better  measure  of  the  ability 
of  the  stream  turbulence  to  promote  transition  than  is  the  total  energy  in  all  frequencies. 

For  Re„  < 2000  the  range  of  unstable  frequencies  (cycles  per  second)  is  given  by 


or 


~ i (0.005  t 0.003) 


(3-1) 


f * (0.005  ± 0.003) 


(3-2) 


The  frequency  which  Is  most  powerfully  amplified  decreases  with  increasing  Re0  so  that  the 
lower  half  of  the  above  frequency  range  may  be  expected  to  be  most  influential  to  transition. 

The  Spangler  and  Wells  experiments  were  carried  out  at  fixed  unit  Reynolds  numbers  in 
air.  For  their  conditions,  the  quoted  unstable  frequency  range  was  20  to  80  cps  which  is  in 
essential  agreement  with  the  above  expression.  Part  of  their  experiment  involved  the  intro- 
duction of  acoustic  disturbances  at  various  fundamental  frequencies  (these  data  are  not  included 
in  Figure  3-3).  Only  those  sound  spectra  with  significant  energy  below  50  cps  were  effective 
in  promoting  transition.  The  Spangler  and  Wells  data  shown  in  Figure  3-3  are  for  grid  turbu- 
lence which  had  relatively  broad  spectra  at  low  frequencies. 
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The  spectra  of  disturbances  measured  by  Schubaur  and  Skramstad  included  significant 
acoustic  energy  at  frequencies  above  50  cps.  This  part  of  the  energy  was  probably  inefficient 
in  promoting  transition,  which  would  tend  tu  explain  the  high  value  of  Re^  T for  I = 0.001 
(in  comparison  to  Spangler  and  Wells,  Reference  3-10). 

Two  simple  alternative  interpretations  of  these  data  are  shown  in  Figure  3-3.  The 
earlier  data  are  reasonably  well  correlated  at  low  disturbance  levels  by 

ReO  J - 64  I'0-"2  (3-3) 

while  the  data  of  Spangler  and  Wells  and  the  blunt  body  data  of  Dunlap  and  Kuethe  are  better 
represented  by 

Re0  T = 6.0  I*0'7  (3-4) 

The  form  of  these  expressions,  that  is,  negative  powers  of  I,  is  chosen  for  simplicity.  These 
expressions  obviously  are  invalid  for  large  1.  There  is  no  way  to  rationalize  the  difference 
between  these  data  trends  in  terms  of  the  single  parameter  1 and  further  investigation  of  the 
role  of  energy  spectra  is  beyond  the  scope  of  the  present  work.  Some  preference  is  felt  for 
the  work  of  Spangler  and  Wells  for  nosetip  application,  especially  since  it  is  in  essentia) 
agreement  with  the  smooth  blunt  body  results  of  Reference  3-7. 

3.2.7  Summary  of  Pertinent  Experimental  Data 

Most  of  the  smooth  wall  transition  data  which  have  been  presented  in  this  section  suffer 
from  a lack  of  accurate  definition  of  stream  disturbances.  The  primary  exceptions  are  the 
data  of  Dunlap  and  Kuet‘,e  and  the  smooth  wall  flat  plate  data.  Information  from  flights  of 
blunt  bodies  in  the  atmosphere  is,  however,  useful  because: 

• It  demonstrates  that  laminar  flow  on  blunt  bodies  can  be  extend-''  very  large 
momentum  thickness  Reynolds  numbers. 

• It  indicates  that  for  smooth  surfaces,  transition  in  the  atmosphere  is  sensitive 
to  minor  surface  irregularities  and,  ■ 'arlations  in  atmospheric  conditions. 

The  latter  point  is  an  indication  that  if  large  uncertainty  in  transition  altitude  below,  say 
40,000  feet,  is  detrimental  to  mission  objectives,  smooth  nosetip  surfaces  should  be  avoided. 
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3.3  TRANSITION  MODELING 


For  smoother  graphite  nosetips,  transition  will  proaably  be  influenced  both  by  surface 
roughness  and  by  atmospheric  disturbances.  It  is  therefore  most  convenient  to  have  a single 
transition  criterion  which  combines  both  effects.  Since  the  influence  of  roughness  is  more 
completely  understood,  the  rough  wall  transition  criterion  is  used  as  a ba  is  for  the  general- 
ized transition  prediction  technique. 


3.3.1  PANT  Rough  Wall  Transition  Criteria 

This  nosetip  transition  prediction  method  was  developed  as  a correlation  of  extensive 
rough  blunt  calorimeter  data  taken  in  a hypersonic  wind  tunnel  at  a free  stream  Mach  number 
of  5.  Model  size,  geometry  and  surface  roughness  height  were  varied  over  quite  wide  ranges. 
This  model  generality,  together  with  variable  tunnel  stagnation  pressure  and  stagnation  tem- 
perature, allowed  simulation  of  the  dimensionless  parameters  Re^,  k/b  and  Tw/Te  for  graphite 
nosetips  during  reentry.  The  data  and  their  correlation  are  described  in  detail  in  Refer- 
ence 3-1. 

The  correlation  is  stated  in  terms  of  the  value  of  an  amplification  parameter.  Re,  at 
transition  as  a function  of  a disturbance  pa-ameter,  k/(^u).  For  the  PANT  data,  j < is  the 
wall  cooling  ratio,  T /T  . The  transition  criteria  are 


Re.  (ft'’ 

>l\  'JiU  / 


( 355  at  = 1 , onset  / 
I 215,  location  ' 


(3-5) 


The  distribution  of  this  transition  parameter  around  a blunt  body  has  its  maximum  near  the 
sonic  point  so  that  the  location  criterion  does  not  predict  transition  in  the  supersonic  region 
between  the  sonic  point  and  the  cone.  Indeed,  in  an  environment  of  increasing  severity,  rough- 
ness induced  transition  "flashes"  forward  from  the  cone  to  somewhere  ahead  of  the  sonic  puirit. 
In  the  PANT  data,  transition  was  not  observed  to  occur  as  far  downstream  as  the  sonic  point. 
This  observation  led  to  the  introduction  of  the  onset  criterion.  The  transition  parameter 
must  exceed  255  at  the  sonic  point  before  the  transition  location  is  determined  as  the  posi- 
tion where  the  parameter  equals  215. 

As  indicated  by  the  transits  n criteria,  the  effect  of  lowering  the  wall  cooling  ratio 
is  to  accentuate  the  influence  of  roughness.  In  Reference  3-1,  arguments  were  presented  which 
indicated  that  the  disturbance  parameter  k//,o  ir.  a measure  of  the  relative  kinetic  energy  at 
the  top  of  the  roughness  elements.  For  non-ideal  gases  it  was  shown,  using  detailed  boundary 
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layer  calculations,  that  > = i:e/f>w  is  more  appropriate  than  ip  = Tw/Te-  The  influence  of  sur- 
face mass  transfer  on  roughness  induced  transition  was  also  interpreted  using  boundary  layer 
solutions.  As  a result,  the  wall  cooling  ratio  parameter  was  generalized  to  read 


, - i!' 
* ' 10 


o * V) 


where  B'  is  the  dimensionless  blowing  parameter 


(3-6) 
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In  the  generalization  of  the  PANT  rough  wall  transition  criteria  to  include  the  smoci.ri 
case  the  momentum  thickness  Reynolds  number  is  retained  as  the  measure  of  the  ability  of  the 
bouno.  ry  layer  to  amplify  disturbances  and  the  disturbance  parameter  is  extended  to  include 
stream  disturbances.  When  both  sources  of  disturbance  are  present  the  total  disturbance  is 
taken  to  be  the  sum  of  the  individual  disturbances,  ear.h  computed  as  if  the  other  were  absent. 
The  additional  information  required  for  the  formulation  of  this  model  is  a smooth  wall  transi- 
tion correlation  for  blunt  bodies.  Since  smooth  wall  blunt  body  data  including  a well  defined 
turbulence  level  is  so  sparse,  the  theoretical  transition  predictions  discussed  in  the  next 
section  are  used  to  help  interpret  the  data. 

3.3.2  Th eoret l r a 1 T r ansitiori  Predictions  for  Smooth  Walls 

Wilcox,  References  3-11  and  3-12,  has  generated  theoretical  transition  predictions  for 
smooth  and  rough  walls  using  a second  order  closure  turbulence  model.  The  essential  features 
of  the  turbulence  are  assumed  to  be  adequately  described  by  the  distribution  of  a specific 
turbulent  energy,  e,  and  a specific  turbulent  dissipation  rate,  it.  The  eddy  viscosity  is 
assumed  to  be  given  by  e * e/a  (which  is  dimensionally  correct)  and  thp  second  order  equations 
for  the  turbulent  fluctuations  are  replaced  by  conservation  equations  for  e and  it.  Several 
closure  constants  appear  in  the  model  and  have  been  previously  evaluated  for  fully  developed 
turbulent  flow.  Two  of  these  were  reevaluated  to  be  Reynolds  number  dependent  in  order  to 
better  predict  smooth  wall  flat  plate  transition. 

These  transition  predictions  for  smooth  flat  plates  agree  well  with  the  trend  of  the 
early  flat  plate  data  discussed  in  Section  3.2.6  (see  rigure  3-4).  For  the  purposes  of  apply- 
ing edge  boundary  conditions,  e = ( 3/2)u ' z . The  value  or  Sz  was  adjusted  to  get  the  best 
agreement  with  these  data. 
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Figure  3-4.  Smooth  wall  transition  predictions  of  Wilcox.  Momentum  thickness  Reynolds 
number  as  a function  of  turbulence  intensity. 


It  seems  correct  to  state  that  the  model  is  tuned  to  agree  with  the  early  smooth  flat 
plate  data.  However,  with  this  tuning  accomplished  and  with  Sifi  held  constant,  the  model  should 
be  c useful  tool  for  extrapolating  to  the  blunt  geometry.  Two  predictions  for  a smooth  hemis- 
phere (Reference  3-12)  are  shown  in  Figure  3-4.  Apparently,  for  these  relatively  large  turbu- 
lence intensities,  the  change  in  configuration  has  little  effect  on  the  transition  Reynolds 
number.  Predictions  for  hemispheres  are  not  yet  available  for  lower  turbulence  intensities. 

It  may  be  that  the  favorable  pressure  gradient  on  the  blunt  geometry  will  be  predicted  to 
delay  transition  at  the  lower  turbulence  levels  (in  comparison  to  the  flat  plate  case).  If 
so,  these  results  would  be  in  disagreement  with  the  data  of  Dunlap  and  Keuthe.  In  view  of 
these  uncertainties,  the  incomplete  results  of  Wilcox  for  smooth  blunt  bodies  are  judged  to 
indicate  that  Re^  ^ as  a function  of  1 may  very  well  be  the  same  on  blunt  bodies  and  flat 
plates. 

3.33  Combined  Disturbance  Mode1 

Since  th‘-  Dunlap  arid  /*<-•♦>, e smooth  wall  blunt  tody  data  favor  the  trend  of  the  Spangler 
arid  Wells  flat  plate  data,  the  srv/>f.h  wall  transition  correlation  is  taken  to  be 

R _T  " 6 l’0'7  (3-8) 

With  this  particular  ralur-  of  the  e/pwot,  this  expression  has  a natural  compatibility  with 
the  rough  wall  transitior  location  correlation,  which  is, 


These  two  expressions  may  now  be  combined  into  a single  transition  criterion  for  use  when  both 
types  of  disturbance-,  are  n.-psent.  The  disturbances  are  taken  to  b'  independent  and  additive 
so  that  a combined  disturbance  parameter  is  defined  as 

/.  ■ ~ + ICO  I (3-10) 

The  transition  criterion  becomes 

He  , j = 210  (3-11  ) 

Predicted  transition  Reynold’,  numbers  are  shown  in  I i gore  j-0  as  a tunc  ion  -A  /!.••  with 
a turbulence  intensity  as  a parameter.  lor  a hninht  rif  ft  I mil  (nprha;:'.  lyjiir.i:  r.f 
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fine  grain  graphites)  on  a nosetip  with  a nose  radius  of  1.5  inches  which  flies  a typical  re- 
entry trajectory,  k/ipo  = 0.3  at  altitudes  below  40,000  feet,  In  that  case,  the  prediction 
shown  in  Figure  3-5  indicates  that  a stream  turbulence  intensity  of  I « 0,001  would  lower 
Re0  T by  about  25  percent.  Taking  a lower  bound  on  roughness  of  0.1  rail,  this  intensity  may 
be  considered  to  be  a rough  lower  limit  on  the  turbulence  intensities  of  interest  to  graphite 
nosetip  transition. 

It  is  not  appropriate  to  use  the  rough  wall  transition  onset  criterion  directly  with  the 
combined  disturbance  model  because  transition  induced  by  stream  turbulence  can  occur  at  any 
location  on  a blunt  body.  In  view  of  the  crude  understanding  of  the  situation,  a rough  approx- 
imation is  suggested  which  gradually  eliminates  the  onset  criterion  as  the  contribution  of  the 
stream  turbulence  increases.  The  onset  criterion  is  rewritten  as 

V0'7  - Conset-  atMe  = 1-°  <3'12> 

where 

i 255(1  - 166  I/&)  | 

Conset  ' ‘"ax  j 215  I (3'13) 

For  166  I/A  > 0.16,  there  is  no  onset  restriction  on  the  application  of  the  location  criterion. 
3.4  FLIGHT  DISTURBANCE  SOURCES 

During  nosetip  reentry  flight,  three  types  of  disturbances,  apart  from  severe  weather, 
seem  most  probable  to  have  an  influence  on  boundary  layer  transition.  These  are  vehicle  vibra- 
tion, clear  air  turbulence  and  cirrus  cloud  ice  crystals.  These  disturbance  sources  have  not 
been  given  in-depth  consideration  in  the  present  study;  however,  some  important  conclusions 
can  be  drawn. 

Data  from  the  flight  tests  of  the  smoothest  calorimeters  reported  in  References  3-3, 

3-4,  and  3-5  are  useful  because  stream  turbulence  levels  can  be  inferred  from  the  combined  dis- 
turbance model  developed  In  Section  3.3.3.  These  inferred  turbulence  levels  are  listed  in 
Table  3-3.  If  these  flights  are  assumed  to  have  taken  place  in  cloudless  skies,  these  in- 
ferred turbulence  intensities  represent  the  influence  of  either  vehicle  vibration  and/or  clear 
air  turbulence  or,  perhaps  for  flight  r-2,  minor  surface  irregularities  as  discussed  in  Sec- 
tion 3.2.2.  The  disturbance  sources  are  discussed  individually  in  the  next  three  subsections. 
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3.4.1  Vehicle  Vibration 


For  the  flight  reported  in  Reference  3-2,  vehicle  vibration  arising  in  the  rocket  engine 
was  judged  to  enhance  the  influence  of  roughness  by  a factor  of  1/3  (see  Section  3.2.1).  It 
is  not  hard  to  visualize  that  a roughness  element  would  be  a more  efficient  boundary  layer  trip 
when  vibrating,  but  vehicle  vibrations  may  not  promote  transition  on  a smoother  surface.  In 
fact,  the  data  from  Reference  3-3  indicate  a very  minor  influence  of  rocket  engine  operation 
on  smooth  nosetip  transition.  Although  the  inferred  disturbance  level  is  increased  by  25  per- 
cent by  the  rocket  engine,  the  Re^  j is  decreased  only  14  percent  to  1 200 . When  compared  to 
very  quiet  wind  tunnels,  this  is  a relatively  high  transition  Reynolds  number. 

The  intensity  of  rocket  engine  vibration  is  expected  to  be  much  greater  than  any  vibra- 
tion experienced  by  a coasting  reentry  vehicle.  Vehicle  vibration,  therefore,  is  judged  to 
have  a very  minor  influence  on  smooth  wall  nosetip  transition. 

3.4.2  Clear  Air  Turbulence 

Clear  air  turbulence  has  been  studied  primarily  in  connection  with  aircraft  operation 
(see  References  3-13  and  3-14,  for  example).  According  to  Reiter  (Reference  3-13),  jet  air- 
craft respond  to  eddies  of  dimensions  of  the  order  of  100  feet.  A reentry  vehicle  passing 
through  this  scale  of  turbulence  at  20,000  ft/sec  will  be  subject  to  disturbances  of  a frequency 
of  200  cps. 

A rough  estimate  can  be  made  of  the  critical  frequency  on  a reentry  body  using  the  ex- 
pression quoted  in  Section  3.2.6  (Equation  (3-2), 

f = 0.005  y (3-14) 

At  an  altitude  of  50,000  feet,  with  Rjj  = 1 inch  and  a velocity  of  20,000  ft/sec,  sonic  point 
values  are  approximately  ufi  = 6,000  ft/sec  and  b » 0.5  rail.  Thus,  f - 7 x 105  or  about  4,000 
times  as  high  as  the  frequencies  expected  from  eddies  of  the  scale  which  affect  jet  aircraft. 

The  eddy  size  corresponding  to  f = 7 x 105  cps  is  about  0.03  foot.  This  is  near  the  scale  at 

which  dissipation  takes  place  and  very  little  energy  is  contained  in  this  range. 

From  this  very  rough  estimate  it  would  seem  that  atmospheric  turbulence  should  not  affect 
nosetip  transition.  Morkovin,  however,  has  expressed  concern  that  there  may  be  enough  energy 
in  the  "tail  of  the  spectrum"  to  trigger  transition  (see  Reference  3-15).  The  inferred  turbu- 
lence levels  listed  in  Table  3-3  would  substantiate  this  concern. 
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Atmospheric  turbulence  is  not  uniformly  distributed  In  any  particular  altitude  layer. 
Around  40,000  to  50,000  feet  altitudes  turbulence  occupies  about  1 percent  of  the  atmosphere 
according  to  a rough  estimate  by  Houbolt  (Reference  3-14).  The  most  complete  characterization 
which  could  be  hoped  for  is  a probability  distribution  with  altitude  of  the  turbulent  spectra. 
Large  uncertainties  in  disturbance  environment  are  indicated. 

Further  problems  are  associated  with  the  analysis  of  the  shock-turbulence  interaction 
and  the  evolution  of  the  stream  turbulence  as  it  passes  around  the  nosetip.  Although  work 
has  been  done  in  this  area  (c.f.,  References  3-16,  3-17,  and  3-18),  it  is  far  from  ready  for 
application. 

It  would  seem  that  the  first  problem  to  resolve  should  be  whether  or  not  clear  air  tur- 
bulence affects  smooth  wall  transition.  Transition  experiments  in  which  smooth  nosetips  are 
flown  through  patches  of  turbulence  of  known  intensity  would  answer  this  question.  If  transi- 
tion is  found  to  depend  on  clear  air  turbulence  intensity,  the  associated  fundamental  uncer- 
tainty in  transition  environment  could  never  be  removed. 

3.4.3  Cirrus  Ice  Clouds 

The  interaction  of  ice  crystals  and  reentry  nosetip  shock  layers  has  been  studied  by 
Finson,  et  al . (Reference  3-19).  Cirrus  clouds  are  made  up  of  ice  crystals  which  are  10  to 
several  hundred  microns  in  diameter.  According  to  the  melt  layer  analysis  of  Finson,  particles 
larger  than  8 microns  survive  the  shock  layer  to  impact  the  surface.  Particles  large-  than  2b 
microns  (1  mil)  impact  the  surface  essentially  unchanged  in  size  with  about  50  percent  of  their 
original  kinetic  energy.  These  particles  leave  impact  craters  several  times  their  size  and 
the  resulting  rough  surface  (k  > 1 mil)  will  cause  transition  at  any  altitude  where  these 
clouds  are  encountered. 

3.5  SUMMARY  AND  CONCLUSIONS 

The  pertinent  smootn  wall  transition  data  base  has  been  reviewed  and  analyzed.  Data 
from  very  smooth  calorimeters  flown  in  the  atmosphere  indicate  that 

t Laminar  boundary  layer  flow  can  occur  on  smooth  blunt  bodies  at  momentum  thickness 
Reynolds  numbers  as  high  as  1400. 

• Transition  on  smooth  surfaces  is  sensitive  to  minor  surface  irregularities  and/or 
variations  in  atmospheric  conditions. 


With  one  shrouded  model  blunt  body  exception,  smooth  wall  transition  experiments  in  which  the 
turbulence  intensity  has  been  documented  have  been  restricted  to  flat  plates.  There  are  dis- 
crepancies among  the  flat  plate  data,  probably  associated  with  the  spectral  distribution  of 
disturbance  energy. 

The  PANT  rough  wall  transition  criteria  have  been  extended  to  include  the  disturbance 
of  stream  turbulence.  The  smooth  surface  transition  model  is  based  upon  the  data  of  Dunlap 
and  Kuethe  (blunt  body)  and  Spangler  and  Wells  (flat  plates)  as  well  as  the  theoretical  analy- 
sis of  Wilcox.  A combined  transition  criterion  based  upon  the  assumption  that  the  disturbances 
are  independent  and  additive  has  been  suggested. 

Data  from  a flight  of  a very  smooth  nosetip,  with  and  without  the  engine  on,  indicate 
that  vehicle  vibration  probably  will  not  affect  smooth  wall  transition.  Clear  air  turbulence 
is  fundamentally  random  in  its  distribution  in  the  atmosphere.  Also,  even  with  well  quantified 
atmospheric  turbulence,  its  effect  on  smooth  wall  transition  is  presently  impossible  to  predict 
and  most  likely  will  remain  so.  Transition  will  almost  certainly  occur  on  nosetips  flown 
through  cirrus  ice  clouds  because  of  surface  cratering. 

Transition  on  smooth  reentry  nosetips  will  be  very  difficult  to  predict  and  may,  in  fact, 
be  fundamentally  uncertain  because  of  random  atmospheric  variations.  If  large  uncertainties 
in  transition  altitude  for  smooth  nosetips  may  cause  thermostructural  failure  or  compromise 
in  some  other  manner  mission  objectives,  then  smooth  nosetips  should  be  avoided.  If  it  is 
important  to  be  able  to  fly  smooth  nosetips,  a study  should  be  conducted  to  assess  the  sensi- 
tivity of  performance  to  large  uncertainties  in  transition  altitude.  If  necessary,  this  could 
then  be  followed  by  further  studies,  perhaps  experimental,  of  the  influence  of  atmospheric 
turbulence  on  smooth  wall  nosetip  transition. 
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f.  1ST  OF  SYMBOLS 


EROS  pressure  gradient  parameter  (c.f.  equation  4-7) 
skin  frUticn  coefficient,  equation  4-9 
heat  transfer  coefficient,  equatic"  4-B 
turbulent  roughness  augmentation  factor 
Doundary  layer  shap;  factor,  aquation  4-3 
stagnation  enthalpy  outside  of  the  boundary  layer 
enthalpy 

mass  diffusion  flux  at  the  wall  of  specie  i caused  by  moleiular  motion  relative 
to  mean  mass  flux 

effective  turbulent  wall  roughness  height 

mass  fret  t ion  of  spei.it  i 

mixing  length 

wall  heat  flu/ 

d'stance  from  symmetry  axis 

i j'hulinf  Reynolds  analo-jv  factor 

Reynold  , ,.r 

wflii-  length  from  ,i<ir|n.ition  point 

ti,i  i , turbu’erl  temper  a lure  I I nutua lion 
treamw  i .o  vr-1ot  i ty 
nor  urn  1 LOmponent  c,i  velocity 
turbulent  velocity  1 1 u-  t o a I ions 
to-  rdin.iLe  iinnim  I t o siiim  ate 

•li'.  taere  from  symmetry  II X i,  t ( . point  whin1  eilqe  1 tr  eanil  ine  trov.r.  lj,e  how  sh'ic  I 


LIST  OF  SYMBOLS  (Concluded) 


8 usual  pressure  gradient  parameter  (see  explanation  of  equation  4-24) 

5 boundary  layer  thickness 

6.|  velocity  defect  thickness,  equation  4-31  or  4-33 

l,  r eddy  viscosities  for  momentum  and  heat 

fn  h 

. viscosity 

c density 

i shear 

u momentum  thickness,  equation  4-2 

Subscripts 

e boundary  layer  edge 

1 inviscid 

u stagnation 


T 


turbot  nut 


4.1 


INTRODUCTION  AND  OVERVIEW 


The  purpose  of  this  task  is  to  assess  the  adequacy  of  the  turbulent  heat  and  mass 
transfer  models  used  in  the  PANT  shape  change  codes  and,  if  necessary,  upgrade  the  prediction 
procedures  to  remove  any  important  Inadequacies.  The  phenomena  of  particular  concern  are 
surface  roughness  and  vortical  layer  effects  on  turbulent  heat  and  mass  transfer  rates.  The 
aspects  of  the  prediction  which  have  been  addressed  include  the  Blasius  type  skin  friction 
law,  use  of  constant  Reynolds'  analogy  factor,  the  boundary  layer  mass  flow  calculation,  lack 
of  consideration  of  vortical  layer  effects  on  the  boundary  layer  profiles,  and  the  use  of 
n:ul  tipi  ica ti ve  augmentation  factors  to  account  for  increased  turbulent  convective  transport 
due  to  roughness.  The  approach  used  is  to  compare  predicted  heat  transfer  distributions  with 
existing  experimental  data  and  exact  theoretical  predictions.  I he  latter  are  needed  because 
of  the  lack  of  an  adequate  data  base  for  the  flight  conditions  of  interest,  namely  high  Mach 
number,  high  Reynolds  number  flows  with  strong  normal  entropy  gradients  and  strong  streamwise 
pressure  gradients. 

The  roles  of  turbulent  transport  in  reentry  vehicle  nosetip  behavior  are  considered  ,n 
this  introductory  section,  along  with  consideration  of  the  factors  which  govern  the  choice  ot 
mathematical  models  to  characterize  the  phenomena.  The  current  methodology  employed  in  the 
PANT  shape  change  codes  is  presented  in  Section  4.2  along  with  some  discussion  of  the  back- 
ground and  significance  of  the  method.  Results  of  the  method  jre  compared  with  experimental 
data  in  Section  4.3,  where  perturbations  to  the  method  are  considered  in  the  interest  of  ob- 
taining improved  correspondence  with  the  axperimental  results.  1 xact  boundary  layer  solutions 
arc  a'  presented  in  an  attempt  to  clarify  the  sources  ot  discrepancies  between  predictions 
di. , u.iid.  It  is  found  in  Section  4.4  that  even  the  exact  solutions  fail  to  predict  vortical 
layer  effect',  due  to  shortcomi ngs  in  current  understanding  of  turbulence  phenomena  under  these 
conditions.  Associated  conclusions  and  recommendations  derived  from  this  study  are  presented 
in  bectluri  4,6. 

4.1.1  The  nole  of  Turbulent  t inn  in  Missile  Nose  tip  Performance  Predict! ons 

The  convection  heat  and  m«-  transfer  are  important  factors  in  determining  the  reces- 
siun  and  thermostructural  rospons'  of  a rei'ritoy  vehicle  nosetip.  Thc  ablation  rate  is  governed 
by  mass  diffusion  a^ro.-S  the  bound  i>  ver  , so  accurate  total  recession  and  shape  change  pre- 
Jictions  depend  on  the  >r  e of  a ri  ad>  turbulent  eor.veo-i vc  transfer  model,  similarly,  tin 

■i  "la  I stress/sir-nn  fields  is  a direct  result  of 


iri-deutl,  thermal  field  which  drives 
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the  convective  heat  transfer  distribution  history.  Therefore,  adequate  predictions  of  nose- 
tip  structural  performance  depend  on  the  technique  used  to  predict  the  turbulent  convective 
heat  and  mass  transfer  distributions. 

Detailed  analyses  of  ground  test  data  and  theoretical  shape  change  predictions  indi- 
cate that  the  total  recession  is  quite  dependent  on  the  history  of  the  shapes  that  develop, 
and  one  can  obtain  substantially  different  total  recession  on  two  models  subject  to  nominally 
identical  environments  depending  on  the  detailed  shape  history  (c.f.  Runs  805  and  815  of 
Series  I Low  Temperature  Ablator  Test  Series,  Reference  4-1).  Presently  it  is  not  clear  what 
local  convective  transfer  rate  prediction  accuracy  is  required  to  be  able  to  predict  total 
recession  and  shape  change  to  within,  say,  20  percent  in  a clear  air  environment.  It  is  pro- 
bable that  the  prediction  should  be  accurate  to  within  20-50  percent  all  of  the  time  and  15- 
20  percent  most  of  the  time. 

4.1.2  Factors  Governing  the  Choice  of  a Prediction  Procedure 

The  procedure  must  be  able  to  adequately  predict  turbulent  transfer  rates  over  a wide 
variety  of  conditions:  from  relatively  smooth  surfaces  to  surfaces  with  scallops  whose  char- 

acteristic length  is  substantially  larger  than  the  boundary  layer  thickness;  for  Mach  numbers 
from  low  supersonic  to  hypersonic;  for  a range  of  shapes  including  slender,  blunt,  concave, 
convex;  for  ablation  rates  typical  of  flight  materials;  etc.  Needless  to  say,  a procedure 
that  can  adequately  account  for  so  many  diverse  effects  must  have  a considerable  degree  of 
flexibility. 

One  of  the  most  important  factors  to  be  considered  in  selecting  a procedure  for  use  in 
a shape  change  code  is  computational  time.  The  PANT  shape  change  codes  are  developed  and  in- 
being  improved  to  serve  as  tools  in  design  and  mission  trade-off  studies.  For  the  codes  to  I <■ 
useful,  one  must  be  able  to  obtain  the  solution  for  a complete  trajectory  in  a few  minutes  on 
a high  speed  digital  computer.  Only  then  will  it  be  possible  to  do  parametric  trade-off  studi*-. 
at  acceptable  cost.  Such  predictions  are  generated  by  marching  through  the  trajectory  in  a 
sequence  of  time  steps,  usually  on  the  order  of  100  to  300.  At  each  time  step  it  is  necessa'y 
to  calculate  the  heat  and  mass  transfer  rate  distributions,  so  an  adequate  procedure  must  Li- 
able to  calculate  these  distributions  r,i  a complete  nosetip  in  a few  seconds,  since  an  e-a<  t 
solution  cf  the  boundary  layer  equations  requires  minutes  of  computer  time  for  the  condit inns 
of  interest,  the  PANT  codes  utilise  a simole  integral  momentum  boundary  layer  solution  a|i|m  i.  r 
with  an  exponential  skin  friction  law  and  boundary  layer  analogy  considerations  to  relate 
local  heat  and  mas  transfer  inefficient-,  to  the  skin  trirtion  coefficient. 


4.2  DISCUSSION  OF  HEAT  AND  MASS  TRANSFER  PREDICTION  MODEL  IN  EROS 
4.2.1  Basic  Nontranspired  Boundary  Layer  Procedure 


Nonblown,  turbulent  convective  transfer  coefficients  in  EROS  are  calculated  utilizing 
a simplified  boundary  layer  integral  momentum  equation  employing  flat  plate  type  friction 
factors  and  Reynolds'  analogy  type  relations.  An  auxiliary  relationship  based  on  a mass  bal- 
ance between  the  boundary  laye.  and  bow  shock  is  used  to  determine  edge  properties  for  a given 
surface  pressure  distribution.  Surface  blowing  effects  are  assumed  to  be  decoupled  from  the 
convective  transfer  calculation  and  are  included  as  multiplicative  factors  after  the  basic 
coefficient  value  is  determined.  The  equations  utilized  are  based  on  the  following  relations: 

• Integral  momentum  equation: 
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where  the  momentum  thickness,  0,  is  defined  by 


and  the  boundary  layer  shape  fat  or,  H,  is  defined  by 


(4-2) 


(4-3) 


The  surface  mass  flux,  (pv)w>  is  taken  as  zero  in  this  decoupled  procedure 

» Crov  ill  "Composite”  skin  friction  correlation  (assumed  to  be  valid  for  all  turbu- 
lent flow  from  a stagnation  point  with  no  surface  mass  transfer,  smooth  wall,  and 
compressible  flow  conditions,  taken  from  Reference  4-3): 
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where  the  "reference  property"  momentum  thickness  Reynolds  numbers  is  defined  as 

ou  n 
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and  the  barred  (or  reference)  properties  are  evaluated  at  the  reference  enthalpy 


h = 0.5  h + 0.3  h + 0.2  H (4-6) 

w e o 

• Reynolds'  analogy  factor 
CH 

Rt  l x--  (=  1.0  according  to  Reynolds'  analogy)  (4-7) 
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2 

where  the  Stanton  number,  Cu,  is 
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and  the  friction  factor,  Cf/2,  is 
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In  practice,  the  mass  transfer  „oef ficier.t 
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is  taken  to  be  equal  to  the  heat  transfer  Stanton  numb-  r by  the  same  arguments  which  justify 
the  application  of  a Reynolds'  analogy  factor  (see  Equations  (4-14)  to  (4-18)  below) 

• Boundary  layer/bow  shock  mass  balance  relationship  (see  subsection  4.2.4) 
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where  the  bow  shock  wave  angle  is  a function  of  yT  which  is  the  distance  from  the  sym- 
metry axis  to  the  point  where  the  'local  edge  streamline  crossed  the  bow  shock.  Mote 
that  y^  can  be  evaluated  exactly  from  the  relation  for  arisymmetric  flow 

Jt  = 2r  ' eue  ( ' * 2/  r<>-v)wd!>  (4-12) 
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Equation  (4-11)  conies  from  neglecting  wall  transpiration  and,  in  part,  from  the  iden- 
tity 

(4  - 6*)  - e (f  - f-)  (4-13) 

whem  the  quantities  and  f-  have  been  evaluated  as  if  the  boundary  layer  were  imcom- 
pressible  with  a l/7th  power  velocity  profile.  Predicted  vortical  layer  effects  on 
heat  transfer  in  later  sections  should  be  considered  in  light  of  these  simplification. 

At  this  juncture  it  is  worthwhile  to  consider  the  boundary  layer  integral  energy 
equation: 
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where  the  energy  thickness,  j,,  is  defined  as 
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The  thermal  boundary  layer  thickness,  4^,  is  in  general  different  than  the  thickness  of  the 
momentum  boundary  layer,  4. 

The  similarity  between  the  momentum  and  energy  integral  equations,  (4-1)  and  (4-14) 
respectively,  Is  the  foundation  for  the  analogy  of  heat  to  momentum  transfer  exemplified  by 
Reynolds  analogy  ( i . e , , the  analogy  come?  from  the  similarity  of  the  boundary  layer  differen- 
tial equations  which  form  the  bases  of  the  integral  equations).  Note  that  for  flow  along  an 
isothermal  flat  plate,  the  bracketed  terms  in  these  equations  are  zero.  If  In  addition  the 
Prandtl  number  is  unity,  the  Reynolds'  analogy  factor  for  this  case  is  exactly  equal  to  unity 
and  the  solutions  of  the  momentum  and  energy  boundary  layer  equations  (integral  and  local)  are 
identical,  assuming  appropriate  definition  of  analogous  terms,  i.e.. 
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Reynolds'  analogy  would  be  expected  to  apply  under  the  conditions  of  unity  Prandtl  number  and 
similarity  in  "acceleration"  terms 


1 + H due  - 1 d(Ho  ~ Hw) 

ue  ds  (Ho  * hwJ  ds 


(4-17) 


Since  in  general  the  similarity  condition  of  Equation  (4-17)  does  not  hold  for  practical  flow 
situations,  Reynolds'  analogy  breaks  down  even  for  unity  Prandtl  number. 

The  composite  wall  shear  law.  Equation  (4-4),  is  based  on  wall  shear  behavior  along  a 
flat  plate.  By  Reynolds'  analogy  (which  is  valid  here  for  Pr  - 1),  it  is  expected  that 
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along  a flat  plate  held  at  constant  temperature.  Accordingly,  it  is  to  be  noted  that  by  re- 
placement of  terms  in  Equation  (4-4)  with  those  presented  in  Equation  (4-18),  a solution  to 
the  energy  integral  equation,  Equation  (4-14),  can  be  obtained  for  non-flat-plate  situations 
Alternatively,  if  H is  taken  as  -1  in  Equation  (4-1)  (as  is  done  in  the  EROS  boundary  layer 
evaluation  methodology)  It  can  be  seen  that  the  solution  that  is  obtained  is,  in  essence,  a 
solution  of  the  energy  equation  for  the  case  of  constant  wall  temperature  and  unity  Prandtl 
number. 

Based  on  these  considerations,  and  anticipating  nonunity  Reynolds  analogy  factors  and 
roughwall  effects  on  heat  and  mass  transfer,  it  is  reasonable  to  formulate  a turbulent  heat 
transfer  relation  (hereafter,  referred  to  as  the  modified  Crowell  composite  heatirg  model) 
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where 
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and  KefJ  is  obtained  by  integration  o!  Equation  (4-1),  with  H = -1.0  and  [pv)w  = 0.  The  de- 
velopment of  Equation  (4-19)  car  be  found  in  Reference  (4-2),  where  the  reasons  for  the  par- 
ticular form  for  bridging  between  the  laminar  and  turbulent  components  of  the  composite 


relation  are  discussed.  The  laminar  Reynolds  analogy  factor,  R^,  and  roughness  augmentation 
factor,  F^,  are  also  discussed  therein  (turbulent  values  will  be  discussed  in  the  next  sut 
sections) . 

Thus,  the  solution  methodology  employed  in  EROS  is  tantamount  to  the  solution  of  the 
energy  integral  relation  for  an  isothermal  wall.  Consequently,  the  "momentum  thickness"  ob- 
tained is  more  likely  the  energy  thickness  which  no  doubt  differs  from  the  momentum  thickness 
by  the  similarity  arguments  presented  above.  Computed  entropy  layer  effects  on  heating  as 
controlled  by  Equation  (4-11)  should  be  judged  in  light  of  this  consideration  (see  Section 
4.3.4). 

4.2.2  Surface  Roughness  Effects  in  Turbulent  F low 

The  increased  turbulent  rates  due  to  surface  roughness  are  predicted  using  a mult'-ii- 
cative  augmentation  factor  which  is  evaluated  using  the  Powars1  correlation,  Reference  4-7. 

^ 1.0  , RKT  r lo  i 

F t = • = I /.oglri( RkTj  + L 10  • RKT  •:  10"/  (4-21) 
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' 3.0  , &KT  > 10"  ' 

The  smooth  wall  Stanton  number,  CH  T>  is  evaluated  from  Equation  (4-19),  taking  FT  = 1.0. 

The  roughness  parameter,  RKT  is: 


RKT  = Rek 
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The  Powar's  correlation  was  developed  for  sand  grain  type  roughness  ranging  from  small 
k k 

(.  1)  to  large  (-- « 1 ) rouqhness  heights.  In  flight,  scallops  form  in  the  turbulent  f’;*  i-ea-.me 

but  the  augmented  convection  does  not  appear  to  be  consistent  with  expects  Ions  considering  scal- 
lop roughness  as  equivalent  to  sand  grain  roughness.  In  thp  PAul  codes,  this  "macroroughness" 
heat  augmentation  is  evaluated  using  an  effect’';.-  "microroughness11  which  is  determined  em- 
pirically for  each  material.  In  ri „,;r  air  environments,  a constant  average  value  of  this 
"turbulent  roughness"  1;  used  after  transition,  and  Powar's  correlation  is  then  used  to 
evaluafo  ; he  instantaneous  augmentation  factor  during  the  calculation.  However,  these  con- 
siderations are  over-ridden  in  weather  environment'  where  computed  crater  depths  due  to 
hydrcmeteor  impacts  are  assumed  to  be  equivalent  to  sand  grain  type  roughness  heights  if  they 
ore  greater  than  the  "turbulent  roughness". 
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4.2.3  Reynolds'  Analogy  Factor  in  Turbulent  Flow 


The  baseline  PANT  codes  utilize  unity  Reynolds'  analogy  factor;  that  is, 

RT  = 1.0  (4-23) 

In  turbulent  heat  transfer  predictions,  the  Reynolds'  analogy  factor  is  often  expressed  as  a 
power  of  the  Prandtl  number,  Ry  « (Pr)a,  where  "a"  may  vary  from  -2/3  to  0.  Such  differences 
in  Reynolds'  analogy  factors  result  in  variations  in  predicted  heat  transfer  rates  of  up  to 
35  percent  for  Pr  = 0.7.  In  assessing  the  most  appropriate  value  for  use  in  the  PANT  codes, 
a review  of  a number  of  exact  numerical  solutions  for  turbulent  boundary  layers  with  zero  pres- 
sure  gradient  and  negligible  entropy  layer  effects  indicated  a value  of  Ky  = (Pr)  1 , Sub- 
sequently, results  of  three  computations  with  streamwise  pressure  gradient  were  analyzed  and 
the  calculated  turbulent  Reynolds'  analogy  factor  exhibited  a strong  dependence  on  pressure 
gradient,  which  is  indicated  in  Figure  4-1.  For  nosetip  applications,  the  pressure  gradient 
parameter,  0,  is  typically  less  than  0.8  or  so,  though  it  can  be  greater  tn  the  vicinity  of 
sharp  corners. 

These  results  indicate  that  the  turbulent  Reynolds'  analogy  factor  may  well  vary  by  as 
much  as  25  to  40  percent  in  nosetip  applications.  Since  predicted  heat  transfer  rate  is 
directly  proportional  to  Reynolds  factor,  this  is  an  important  effect.  To  evaluate 

the  adequacy  of  the  current  procedure,  EROS  code  predictions  were  made*  with  RT  ::  1.0  as  well 
as  with 


i 1.126  , A<0. 15 

Rt  = ' 0.74A'°*22,  0.15  < A ■ 0.32  (4-24) 

/ 

' 0.95  , A=-0. 32 

where  A * (^)  is  a measure  of  the  pressure  gradient.  (The  expression  A is  used  iri  uiate 

wt  2 f ^ue 

of  the  usual  pressure  gradient  parameter,  6 - — where  •;  is  the  Levi-Lees  streamwise  co- 

ue 

ordinate,  since  f,  is  not  available  in  ‘.he  EROS  code.)  The  cutoff  Ry  = 0.95  for  A ■ 0.32  was 
made  after  preliminary  comparisons  witi  PANT  Series  B data  (Reference  4-4)  indicated  that 
lower  values  of  Ry  for  A > 0.32  resulted  in  predicted  heat  transfer  rates  which  are  much  too 
low.  A plot  of  the  pertinent  calculations  and  the  approximate  curvefit  used  for  Ry  in  the 
EROS  calculations  (Equation  (4-24))  i;  shown  in  Figure  4-2. 

Results  of  these  calculations  are  presented  in  Section  4.3. 
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Figure  4-1.  Theoretical  turbulent  Reynolds  analogy  factor  dependence 
upon  pressure  gradient. 
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Figure  4-2.  Curve  fit  for  theoretical  variable 'turbulent 
Reynolds  analogy  factor. 
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Referring  now  to  the  last  paragraph  of  Subsection  4.2.1  above,  it  is  probable  from  a 
theoretical  point  of  view  that  Reynolds  analogy  factors  other  than  Pr"  = 1.126  need  not  be 
applied  to  Equation  (4-19)  since  the  solution  methodology  already  accounts  approximately  for 
nonsimilar  effects  on  the  energy  boundary  layer.  That  is,  if  application  of  Reynolds'  analogy 
is  at  all  valid,  it  should  be  in  reverse  in  order  to  estimate  the  friction  factor  from  the 
computed  Stanton  number.  This  reasoning  provides  an  explanation  for  the  aforementionea  need 
to  limit  the  decreases  in  R-j.  with  increasing  free  stream  acceleration  (Equation  (4-24)). 

4.2.4  Vortical  Layer  Effects  in  Turbulent  Flow 

As  the  boundary  layer  grows  along  the  body  surface  it  ingests  or  swallows  fluid  which 
has  passed  through  weak  as  well  as  strong  portions  of  the  curved  bow  shock.  Since  the  entropy 
behind  the  weak  portion  of  the  shock  is  much  less  than  that  behind  the  strong,  normal  portion, 
an  entropy  gradient  or  "entropy  layer"  exists  behind  the  jhock.  The  term  entropy  swallowing 
refers  to  the  ingestion  by  the  boundary  layer  of  fluid  of  varying  entropy  level  (the  ingestion 
relationship  used  in  EROS  in  Equatio’  (4-11)).  Since  entropy  gradients  can  be  related  to 
vorncity,  this  phenomenon  has  a’so  been  termed  the  "vertical  layer"  effect.  The  impact  of 
“entropy  swallowing"  on  the  turbulent  boundary  layer  is  two-fold.  First,  tne  fluid  which  has 
assec  through  weaker  portions  of  the  shock  will  have  a higher  total  pressure  resulting  in 
higr.er  boundary  layer  edge  momentum  flux  for  a given  local  statu  pressure.  Second,  the  vor- 
ticity  or  -edge  velocity  gradients  will  distort  tne  boundary  layer  profiles  from  their  usual 
constant  entropy  form.  Also  the  edge  vorticity  probably  interacts  with  the  boundary  layer 
turbulent  eddies  resulting  in  altered  boundary  layer  tmnsport  properties.  The  net  effect  of 
these  phenomena  is  to  increase  the  heating  rate  over  that  which  would  be  obtained  if  all  of 
the  flow  has  passed  through  a normal  shock.  The  upper  heating  limit  corresponds  approximately 
to  a cone-flow  oblique  shock  entropy  situation.  The  magnitude  of  this  change  from  normal  to 
oblique  shock  heating  rates  depend  on  many  factors,  including  shock  Lurvature,  free  stream 
Mach  number,  and  boundary  layer  growth.  It  has  been  previously  shown  (Reference  4-5)  that  for 
conditions  typical  of  RV's,  the  effect  can  be  on  the  order  of  a factor  of  two  increase  at 
critical  nosetip  locations. 

The  current  tR05  model  for  predicting  turbulent  heat  transfer  does  not  have  the  flexi- 
bility to  account  for  entropy  layer  effects  on  the  boundary  layer  profiles,  that  is  the  histor- 
ical and  local  effects  of  edge  vorticity  on  the  wake  re.jion  of  the  boundary  layer.  However, 
these  considerations  could  be  incorporated  i-  the  model  through  appropriate  modifications  of 


y 
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the  Stanton  number  relationship,  Equation  (4-19),  based  on  correlation  either  of  experimental 
heat  transfer  data  or  exact  numerical  boundary  layer  results.  An  alternate  approach  to  mod- 
eling the  behavior  was  considered;  at  any  given  streamwise  station  at  which  entropy  layer 
effects  are  important  the  local  edge  entropy  employed  in  the  solution  of  the  momentum  inte- 
gral equation  is  considered  to  be  the  mass  weighted  average  of  the  inviscid  flow  entropy 
which  has  been  entrained  by  the  boundary  layer  to  that  station.  The  effect  of  the  mass  aver- 
aging of  entropy  is  to  reduce  the  heat  flux  predicted  by  the  EROS  procedure.  Comparisons  with 
data  indicate  that  the  mass  averaged  entropy  layer  procedure  results  in  substantially  better 
predictions  at  high  Mach  number  conditions  (see  Subsection  4.3.3). 

4.3  COMPARISON  OF  EROS  PREDICTIONS  WITH  EXPERIMENT  AND  "EXACT'1  SOLUTIONS 

In  this  subsection,  EROS  predictions  of  heat  transfer  coefficients  are  compared  to 
experimental  data  for  a variety  of  body  shapes  and  free  stream  conditions,  including  wind 
tunnel  and  ballistic  range  environments.  The  ballistic  range  heat  transfer  coefficient  data 
are  obtained  by  inference  of  heat  flux  from  measurement  of  a melt  progression  line  on  a metal 
calorimeter  model*.  Although  the  heat  transfer  data  are  indirect  and  subject  to  uncertainty, 
the  data  are  important  because  the  environments  match  actual  flight  conditions  more  closely 
than  do  wind  tunnel  facilities. 

4.3.1  Wind  Tunnel  Data  Comparisons,  M = 5 Smooth  Wall 

The  smooth  wall  heat  transfer  coefficient  data  used  in  this  comparison  were  obtained 
during  the  PANT  Series  8 tests  in  the  Naval  Surface  Weapons  Center  (NSWC)  wind  tunnel  number 
8 at  M = 5.  A full  description  of  the  calorimeters  and  test  conditions  can  be  found  in 
References  4-4  and  4-6.  In  order  to  isolate  any  inadequacies  in  the  EROS  boundary  layer  pre- 
diction technique,  a combination  of  results  from  experimental  measurements  and  exact  inviscid 
calculations  (using  the  RAZIBB  code.  Reference  4-7)  were  utilized  to  establish  the  correct 
surface  pressure  and  shock  shape  data  for  the  boundary  layer  calculations.  Also,  experimen- 
tally determined  wall  temperature  distributions  were  input  into  EROS. 

Four  predictions  were  made  using  EROS  for  each  of  the  three  Series  B models: 


See  Reference  4-5  for  a discussion  of  the  data  reduction  technique. 


4-15 


Case  1.  Baseline  EROS  prediction  which  Includes: 

• Modified  Crowell  composition  heating  law  (Equation  (4-19)) 

• Unity  turbulent  Reynolds*  analogy  factor,  = 1.0 

• Baseline  entropy  layer  swallowing  calculation  (c.f.  Section  4.2) 

Case  2.  Same  as  1 except  with  streamwise  pressure  gradient  dependent  Reynolds'  analogy 
factor.  Equation  (4-24). 

Case  3.  Same  as  1 except  with  thr  boundary  layer  edge  conditions  being  calculated 
using  the  mass  average  model  described  in  Subsection  4.2.4. 

Case  4.  Baseline  EROS  modified  with  2 and  3 

In  addition,  selected  exact  turbulent  boundary  layer  solutions  were  generated  using 
the  BLIMP  code*  with  three  different  models  for  the  turbulent  eddy  viscosity:  Kendall  (Ref- 

erence 4-8),  Bushnell /Beckwith  (Reference  4-9),  and  Cebeci/Sini th  (Reference  4-10).  Experi- 
mental wall  temperatures  were  also  input  to  BLIMP. 

Comparisons  of  predictions  with  the  data  for  the  simple  biconic,  convex  biconic,  and 
triconic,  models  are  shown  in  Figure  4-3.  4-4  and  4-5  respectively.  All  BLIMP  predictions  in 
Figures  4-3  to  4-5  are  with  Kendall's  turbulent  transport  properties.  Important  observations 
from  these  comparisons  are  presented  in  the  following  paragraphs. 

Simple  Biconic  (Figure  4-3) 

• All  predictions  underpredict  the  heat  transfer  coefficient 

• Best  agreement  with  data  among  the  EROS  predictions  is  obtained  with  the  modified 
Crowell  composite  heating  model  and  variable  Reynolds’  analogy  factor  (Case  2). 

Note  that  the  pressure  is  essentially  constant  on  the  45”  face,  so  the  pressure 
gradient  dependent  Reynolds'  analogy  factor  for  this  case  is  constant  at  R,  = 1.126 
(i.e.,  Pr  'A).  The  enmposi t.e/heating  model  gives  good  agreement  in  shape  of  the 
heat  tran'fjr  coefficient  distribution,  but  the  level  is  about  7 percent  low. 


Boundary  Layer  Integral  Matrix  Procedure,  a code  for  obtaining  the  exact  numerical  solution 
of  the  comolete  boundary  layer  equations  for  laminar  or  turbulent  flow,  including  the  effects 
of  edge  vorticity,  mass  addition,  transverse  curvature,  etc.  (c.f.,  Reference  4-13). 
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• Implementing  the  approximate  mass-averaged  vortical  layer  effects  reduces  the 
nominal  heat  transfer  coefficient  by  about  15  percent,  leaving  a disagreement  be- 
tween theory  and  data  of  about  22  percent  with  R-j.  = 1.126  and  30  percent  with  R-j.  * 

1 .0 

• On  the  aft  cone  sidewall,  best  agreement  is  obtained  with  variable  Ry.  Note  that 
there  is  no  noticeable  effect  of  the  approximate  mass  averaged  vorticity  effects 
model  in  this  region. 

• The  BLIMP  prediction  agrees  well  in  shape  and  is,  on  the  average,  about  5 percent 
low. 

Before  looking  at  additional  comparisons  with  data,  a few  cornnents  about  the  BLIMP 
predictions  are  in  order.  At  this  point  it  is  well  to  note  that  typical  eddy  viscosity  models 
depend  on  terms  like  the  boundary  layer  thickness  or  mean  flow  normal  velocity  gradients  which 
are  presumed  to  vanish  at  the  boundary  layer  edge.  For  boundary  layers  in  the  presence  of 
inviscid  entropy  gradients,  there  is  no  currently  available  unambiguous  definition  of  boun- 
dary layer  thickness,  and  the  mean  velocity  gradient  does  not  vanish.  Recently,  the  BLIMP 
code  was  modified  to  include  the  option  of  utilizing  two  other  turhulent  eddy  viscosity  models 
Bushnell/Beckwith  and  Ceheci/Smi th.  Predictions  of  the  simple  biconic  data  (Figure  4-3)  with 
these  two  models  are  compared  with  the  predictions  using  Kendall's  model  in  Figure  4-6.  Ig- 
noring for  the  moment  the  curve  on  this  figure  for  the  Kendall  model  with  the  velocity  con- 
straint, it  is  evident  that  all  turbulent  models  do  a good  job  of  predicting  the  shape  of  the 
heat  flux  distribution,  with  the  Kendall  model  doing  the  best  job  of  predicting  the  levels  as 
well.  However,  the  success  of  the  shapes  is  attributable  both  to  the  apparent  approximate 
validity  of  all  of  the  turbulent  models,  and  to  a computational  manipulation  to  lend  stability 
to  the  solutions.  Without  this  manipulation,  the  predicted  average  levels  are  reasonable,  but 
the  shapes  are  poor,  as  noted  in  the  curve  for  the  Kendall  model  with  the  velocity  constraint. 
The  difficulty  lies  in  the  finite  nature  of  the  velocity  gradient  at  the  edge  of  the  boundary 
layer  which  causes  the  values  of  eddy  viscosities  in  the  outer  portiors  of  the  boundary  layer 
as  presently  modeled  to  be  sensitive  to  the  choice  of  location  of  the  boundary  layer  edge, 
at  least  for  the  Kendall  and  Cebeci  models.  The  trick  that  was  used  to  circumvent  this  dif- 
ficulty was  to  base  the  definition  (by  the  code)  of  the  edge  of  the  boundary  layer  on  the 
behavior  in  the  boundary  layer  on  the  stagnation  enthalpy  profile  that  has  no  gradient  at  the 
edge.  A more  complete  discussion  of  this  problem  is  given  in  Subsection  4.4,  along  with  a 


recommended  solution.  All  BLIMP  results  other  than  the  one  curve  shown  in  Figure  4-6  have  been 
based  on  the  enthalpy  profile  constraint  as  a means  to  dampen  the  erratic  nature  of  solutions 
with  current  eddy  viscosity  models. 


Convex  B iconic  (Fiqrre  4-4) 

• All  of  the  EROS  predictions  are  adequate  on  the  torebody,  except  that 

- near  the  nose,  with  variable  R^.,  the  heat  transfer  coefficient  is  overpredicted 
by  about  13  percent 

- all  of  the  procedures  overpredict  heat  transfer  near  the  shoulder.  This  is  pro- 
bably a result  of  inaccuracies  in  the  local  pressure  as  a result  of  uncertainty 
in  the  local  body  slope, 

• Base  EROS  is  best  on  the  aft  cone  (about  B percent  high)  with  the  variable  R-j. 

resulting  in  predictions  about  17  percent  high. 

• BLIMP  substantially  overpredicts  the  heating  on  the  forebody 


Triconic  (Figure  4-5) 

• Overall,  the  base  EROS  prediction  is  best  on  the  forebody  and  is  adequate  on  the 
aft  cone,  though  improvement  is  obtained  by  inclusion  of  the  approximate  mass 
averaged  entropy  layer  effects  on  the  aft  cone. 

• On  the  forecone,  inclusion  of  the  mass  averaged  entropy  layer  effects  results  in 
the  predicted  maximum  heat  transfer  coefficient  being  low  by  about  25  percent. 

• Variable  Ry  results  in  predictions  which  are  about  9 percent  high  at  peak  heating 
and  13  percent  high  on  the  aft  cone. 

t Inclusion  of  variable  R^.  and  mas;  averaged  entropy  layer  effect;  yields  a peak 
heating  prediction  which  is  about  20  percent  low  at  peak  heating. 

• The  BLIMP  prediction  results  in  very  good  agreement  with  the  data  except  in  the 
vicinity  of  the  shoulder. 

4.3.2  Rough  Wall  Wind  Tunnel  Data  Comparison,  M = 5 


i 
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The  four  EROS  predictions  are  compared  with  the  PANT  Series  0 3.5  mil  roughened  60°  < 
biconic  data  in  Figure  4-7.  The  predictions  were  made  with  Powars1  roughness  augmentation  ’ 
correlation,  Equations  (4-21)  and  (4-22).  Also  shown  are  smooth  wall  predictions  using  the  i 
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BLIMP  code  and  the  base  EROS  procedures.  None  of  the  predictions  shows  good  agreement  with 
the  trend  of  the  data,  though  the  base  EROS  plus  approximate  mass  averaged  entropy  layer  ef- 
fects prediction  agrees  well  with  the  average  level  of  heat  transfer  coefficient  on  the  face 
of  the  model.  Also  note  that  inclusion  of  the  pressure  gradient  dependent  variable  turbulent 
Reynolds'  analogy  factor  results  in  very  poor  agreement  in  the  shape  of  the  heat  transfer 
coefficient  distribution.  These  predictions  are  based  cn  the  experimental  shock  shape  and 
Dahm/Love  pressure  distribution  correlation  (which  is  in  the  PANT  codes).  Since  this  con- 
figuration is  similar  to  those  for  which  it  has  been  demonstrated  that  the  Dahm/Love  correlation 
yields  quite  good  results  (Reference  4-11),  the  poor  agreement  in  shape  is  probably  not  a 
result  of  utilizing  the  correlation  for  pressure  distribution.  Rather,  it  probably  indicates 
that  the  procedure  for  correlating  Ry  versus  pressure  gradient  is  inadequate  (c.f.  previous 
discussion  at  the  end  of  Subsection  4.2.3). 

Finally,  it  is  useful  to  consider  the  smooth  wall  predictions.  The  BLIMP  prediction 
shows  good  qualitative  agreement  with  the  data  on  the  latter  half  of  the  forecone  if  one  pre- 
sumes that  the  increased  heating  due  to  roughness  is  approximately  a multiplicative  constant 
over  the  smooth  wall  value.  The  BLIMP  and  base  EROS  predictions  agree  well  qualitatively, 
with  the  BLIMP  prediction  being  about  12  percent  higher  than  the  EROS  prediction,  a difference 
which  is  fairly  common  in  predictions  of  turbulent  heat  transfer.  The  trend  of  departure  of 
the  data  from  the  BLIMP  predictions  is  indicative  of  a greater  entrainment  by  the  boundary 
layer  of  fluid  coming  through  the  oblique  portion  of  the  shock.  A greater  entrainment  would 
be  expected  considering  the  additional  shear  due  to  roughness. 

4.3.3  Ballistic  Range  Pat..  Comparisons,  M - 16 

Comparisons  between  predicted  and  inferred  heat  transfer  rates  on  melting  calorimeter 
models  in  the  AEDC  ballistic  range  arc  given  in  Eiguias  4-8  3nd  4-9.  Both  models  are  flat 
faced  45“  biconics.  Best  estimates  of  the  Shock  shape  and  surface  pressure  distributions 
were  obtained  by  synthesizing  experimental  shadowgraphs  and  exact  solutions  (c.f.,  Sandhu  and 
Laub,  Reference  4-5,  for  details).  Four  different  predictions  are  compared  with  the  inferred 
heat  transfer  coefficient: 

1.  SAANT*  predictions  with  the  old  skin  friction  law  and  Ry  = 1.0  (see  Reference  4-12)  ■ 

2.  Base  EROS  with  modified  Crowell  composite  heating  law  and  Ry  = 1.126 

★ 

Steady  State  Analysis  of  Ablating  Nosetips,  an  early  version  of  the  EROS  code,  c.f.  P.cfcrcncc 
4-18. 


4-24 


% feMA/T*) 


omparisor.  of  analytically  predicted  heating  levels  to  those  inferred  from 
allistic  range  data.  Shot  No.  -620. 
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Figure  4-9.  Comparison  of  analytically  predicted  heating  levels  to  those  inferred  from 
ballistic  range  dcta.  Shot  No.  3618. 


3.  Same  as  2 except  approximate  mass  averaged  entropy  layer  effects 

4.  B„IMP  predictions  using  Kendall's  turbulent  ecdy  viscosity 

Examination  of  the  results  shows  that  the  SAANT  and  base  EROS  procedures  with  Rj  = 1.126 
grossly  overpredict  the  heat  transfer  coefficient,  the  latter  by  more  than  100  percent  for 
shot  3618,  However,  inclusion  of  the  mass  aver'.ged  entropy  layer  model  results  in  much  better 
agreement  with  the  data,  the  predictions  neing  within  10  to  25  percent  of  the  data  and  agreeing 
very  well  with  the  shape  of  the  heat  transfer  coefficient  data.  The  BLIMP  prediction  agrees 
moderately  well  with  the  data  for  these  two  shots. 

4.3.4  Discussion  of  Comparisons 

The  above  comparisons  must  be  judged  while  keeping  in  mind  the  potential  inaccuracies 
of  the  heat  transfer,  shock,  and  pressure  distribution  data;  the  inaccu»"’~ies  in  evaluating 
the  theoretical  models  (i.e,,  numerical  inaccuracies);  and  weaknesses  n the  theoretical 
models  themselves  due  to  poor  approximations  of  the  physics.  It  is  elieved  reasonable  to 
assume  that  at  least  for  the  NSWC  experiments  the  heat  transfer  data,  shock  data,  and  pres- 
sure distribution  data  may  be  considered  accurate  based  on  the  exactitude  of  the  experiments 
and  past  data  evaluations  and  the  availability  of  good  shock  and  pressure  data  from  both  ex- 
periment and  exact  inviscid  flow  calculations.  Thus,  in  what  follows,  departures  of  predic- 
tions from  experimental  data  are  taken  to  be  attributable  to  a combination  of  numerical  and 
theoretical  modeling  shortcomings. 

In  the  case  of  results  of  numerical  boundary  layer  solutions  such  as  obtained  from 
BLIMP,  a certain  degree  of  error  will  exist  for  turbulent  flow  situations  because  of  inade- 
quate characterization  of  turbulence  phenomena  which  ultimately  manifest  themselves  in  the 
forms  of  eddy  viscosity,  eddy  conductivity,  and  eddy  mass  diffusivity.  These  are  influenced 
by  all  manner  of  things  including  wake  intermittency,  wall  roughness,  transitional  phenomena, 
and  so  on.  However,  given  good  numerics,  it  is  comforting  to  know  that  one  need  not  be  con- 
cerned about  effects  of  such  things  as  compresibility,  acceleration,  wall  streamwise  temper- 
ature gradients,  reference  properties,  Reynolds'  analogy  factors,  continuity  considerations, 
and  so  on,  ad-infinitem,  that  must  be  considered  with  integral  approaches.  The  necessary 
considerations  to  implement  successfully  the  integral  techniques  are  determined  naturally 
within  the  "exact"  numerical-boundary  layer  solution  procedures  where  the  results  in  theory 
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depart  from  experiment  only  by  the  shortcomings  of  the  physical  modeling  and  the  limitations 
of  classical  boundary  layer  tneory.* 

Two  important  Known  weaknesses  in  BLIMP  modeling  (in  addition  to  transitional  modeling 
shortcomings)  include  the  effects  of  free  stream  vorticity  on  wake  eddy  diffusivities  and  in- 
termittency.  However,  these  latter  models  may  be  adequate  with  relatively  minor  extensions 
in  their  present  forms.  Some  known  and  previously  mentioned  numerical  difficulties  which  are 
discussed  in  the  next  subsection  have  limited  the  assessment  of  current  model  appl icaoi 1 i ty 
to  vertical  flows. 

The  factors  governing  the  choice  of  a prediction  procedure  have  been  outlined  in  Sub- 
section 4.1.2.  The  computational  time  constraints  dictate  the  desirability  of  integral  boun- 
dary layer  solution  techniques.  It  is  worth  emphasizing  that  the  integral  boundary  layer 
equations  are  as  exact  as  the  nonintegral  equations,  and  solutions  can  be  forced  to  conform, 
given  latitude  for  adjustments  of  the  laws  and  “constants"  Incorporated  into  the  integral 
equations. t 

It  can  be  concluded  that  the  present  baseline  EROS  procedure  does  very  well,  considering 
all  of  the  phenomena  that  occur  for  the  test  conditions  selected  herein  for  comparison  pur- 
poses. However,  the  most  significant  failure  of  the  overall  technique  is  unfortunately  for 
the  test  ;onditions  most  applicable  to  the  flight  environment,  i.e.,  the  comparisons  with  the 
ballistic  range  data  of  Figures  4-3  and  4-9.  The  most  significant  reason  for  this  failure 
can  be  derived  by  expanding  on  the  discussion  at  the  end  of  Subsection  4.2.1.  Specifically, 
it  is  believed  that  the  "entropy  swallowing"  relation  ( Eqi  at i on  (4-11))  is  both:  (1)  an  in- 

adequate representation  of  the  exact  relation  (Equation  (4-12));  and  (2)  is  improperly  eval- 
uated in  the  current  procedure. 

For  the  conditions  considered  herein  (no  transpiration)  the  assei  .ion  of  the  inadequacy 

of  the  relation  stems  from  the  evaluations  employed  for  the  boundary  layer  shape  factors.  It 

is  known  that  for  a boundary  layer  with  a power  law  velocity  profile,  the  entropy  swallowing 
6-6* 

shape  factor,  — , increases  with  increasing  edge  Mach  number  and  wall  temperature.  Thus, 

★ 

For  example,  it  can  be  argued  that  wall  roughness  (vis-a-vis  waviness)  is  such  as  to  violate 
the  applicability  of  the  boundary  layer  equations. 

.1. 

From  a practical  standpoint,  the  integral  relations  are  in  some  cases  superior  to  working  with 
the  detailed  boundary  layer  relations.  For  example,  experimental  wall  roughness  effects  can 
be  easily  incorporated  into  an  integral  solution  without  conceptual  difficulty,  unlike  the 
problem  of  detailed  modeling  in  the  non-integral  relations. 
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the  relation  is  at  least  In  error  by  lack  of  consideration  of  these  effects  (the  errors  broaden 
for  the  flight  case  with  ablation).  It  is  probable  that  for  the  conditions  of  these  experi- 
ments, the  entropy  swallowing  shape  factor  is  less  than  that  implicit  in  Equation  (4-11) 
because  of  the  low  wall  to  edge  temperature  ratios,  implying  that  the  evaluated  magnitude  of 
yT  is  too  large  due  to  the  implicit  approximations. 

Concerning  the  assertion  of  improper  evaluation  of  Equation  (4-11),  the  values  of  Re(i 
which  are  inserted  into  the  relation  are  probably  too  large  (also  making  y^  too  large)  for 
two  reasons: 

1.  The  composite  wall  shear  law  is  such  as  to  increase  Re.,  at  a given  body  location 
from  its  actual  value  for  all  cases  where  the  flow  is  not  turbulent  over  the  entire 
run  length.  It  can  be  seen  from  the  form  of  Equation  (4-11)  that  momentum  thick- 
ness increases  (but  not  linearly)  as  the  average  wall  shear  increases  upstream  of 

a given  point.  Smooth  wail  laminar  and  transitional  wall  shears  are  less  than 
those  represented  by  the  composite  law,  and  these  lower  values  do  in  fact  exist 
upstream  of  body  points  with  fully  developed  turbulent  boundary  layers. 

2.  More  important,  the  "momentum  thickness"  derived  from  the  procedure  is  more  properly 
identified  as  the  energy  thickness  (which  does  not  in  general  scale  with  the  boun- 
dary layer  thickness  i.e.,  the  entropy  swallowing  shape  factor  Is  in  error).  It 

is  expected  that  the  momentum  thickness  for  geometr._s  and  conditions  of  interest 
will  be  found  to  be  less  than  the  energy  thickness.* 

Thus,  all  theoretical  considerations  point  to  evaluated  values  of  yT  which  are  too  large. 
Consequently,  the  entropy  layer  is  swallowed  too  quickly  in  all  cases  considered.  The  suggested 
effect  of  the  composite  shear  law  is  demonstrated  dramatically  in  comparing  the  FROS  results 
with  the  SAANT  results  In  Figures  4-8  and  4-9  (i.n.,  peak  heating  is  shifted  forward  on  the 
body  and  to  a higher  level  using  the  composite  law).  In  addition,  the  effort  to  improve  com- 
parisons by  averaging  the  "swallowed"  entropy  is  equivalent  to  demonstration  of  the  need  to 


For  example,  integrating  simultaneously  the  momentum  and  energy  integral  equations  for  an  iso- 
thermal wall,  and  employing  the  Chllton-Colburn  analogy 


£ Or2/  «'  2 (ranges  from  about  0.5 

4>  '3  13  + n to  0.7  for  Pr  =■  0.7) 


The  exact  BLIMP  solution  at  the  stagnation  point  of  the  simple  hiconic  at  Figurp  4-.T  and  4-6 
yields  a value  of  e/4  = 0.68. 
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lower  yy.  And  finally,  BLIMP  evaluates  Equation  (4-12)  exactly  (although  in  terms  of  stream 
function  rather  than  boundary  layer  thickness  parameters)  and  BLIMP  does  a good  job  of  pre- 
dicting the  ballistic  range  data  with  no  need  for  remodeling  either  the  physics  or  the  numer- 
ical procedures. 

In  order  to  improve  the  accuracy  of  the  integral  solution  approach  while  retaining  its 
simplicity,  the  procedure  should  be  upgraded  to: 

1.  Solve  the  full  integral  momentum  equation  ( i . e . , include  H)  to  at  least  provide  an 
improved  entropy  swallowing  evaluation  and,  at  the  upper  end  of  sophistication 
possibly  provide  velocity  profile  information  to  be  used  in  the  energy  equation 
integral  parameters. 

2.  Solve  the  energy  integral  equation  to  at  least  provide  good  heat  transfer  information 
(and  mass  transfer,  by  analogy)*  given  improved  edge  condi tions(though  improved 
entropy  swallowing)  and  at  the  upper  end  of  sophication  possibly  provide  density 
profile  information  for  evaluation  of  momentum  integral  functions  (perhaps  mini- 
mizing the  need  for  reference  properties). 

3.  Modify  the  required  heat  and  momentum  transfer  lows  to  account  for  transpiration, 
roughness,  acceleration,  transition,  and  vorticity. 

As  the  formulation  currently  stands,  it  considers  neither  of  the  "acceleration"  terms 
presented  in  Equation  (4-17).  Both  of  them  are  important,  although  the  streamwise  wall  tem- 
perature gradient  effects  are  certainly  the  lesser  of  the  two  (no  doubt  the  reason  why  the 
present  method  does  so  well  when  entropy  layer  effects  are  minimal). 

It  can  be  argued  that  if  the  entropy  swallowing  calculation  is  improved,  EROS  type 
predictions  will  yield  lower  heat  transfer  coefficients  for  the  NSWC  data  and  thus  render  com- 
parisons relatively  unfavorable.  Three  considerations  suggest  otherwise: 

1.  Entropy  layer  effects  are  not  nearly  so  significant  at  M - 5.  For  example,  at  a 
45”  shock  angle,  the  ratio  of  total  pressure  at  the  stagnation  point  to  that  behind 
the  oblique  shock  decreases  from  about  0.3  at  M - 5 to  about  0.2  at  M = 20.  Thus, 
for  the  same  shock  shape,  boundary  layer  edge  velocities  and  densities  do  not  in- 
crease as  muc.i  at  M = 5 for  a completely  swallowed  vortical  layer  as  they  do  at 

★ 

The  mass  boundary  layer  is  more  analogous  to  the  energy  boundary  layer  than  energy  is  to  the 
momentum  boundary  layer. 
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M * 20.  Secondly,  shocks  become  more  blunt  for  a given  body  shape  as  the  Mach 
number  is  reduced,  moving  entropy  swallowing  further  downstream  on  the  body, 

2.  The  Reynolds'  analogy  factor  that  should  be  employed  in  the  solution  of  the  energy 
integral  equation  probably  ought  to  be  greater  than  1.0  (like  P “V3)  and  probably 
invariant  for  accelerating  flows.  This  suggestion  is  based  both  on  theoretical 
concepts  and  on  the  observation  that  the  variation  of  Rj  employed  herein  invari- 
ably did  a disservice  to  the  streamwise  variation  of  predicted  heat  transfer  coef- 
ficient in  every  case  where  it  departed  from  its  maximum  value  (=  Pr~ V» ) near  the 
nose. 

3.  Heat  transfer  coefficient  results  will  change  somewhat  when  wall  streamwise  temper- 
ature gradients  are  considered. 

Two  approaches  to  the  improvement  of  the  heat  and  momentum  transfer  laws  might  be  con- 
sidered: 

1.  Modify  the  wall  shear  and  wall  heat  flux  laws  directly  to  account  for  all  the  ef- 
fects considered  above.  It  is  suggested  that  even  though  the  composite  heating 
law  has  yielded  some  improvement  in  predicted  heat  flux  distributions,  the  law  is 
too  cumbersome  when  transpiration,  roughness,  acceleration,  transition,  and  vor- 
ticity  are  considered  simultaneously.  Furthermore,  the  real  physics  of  turbulent 
starting  length  and  transitional  effects  are  violated  with  the  corposite  model. 

2.  A more  sophisticated  approach  is  to  extend  the  two-parameter  methods  for  predicting 

skin  friction  (e.g..  Reference  4-16)  in  which  velocity  profiles  are  described  in 
terms  of  a "law  of  the  wall"  and  a "law  of  the  wake"  with  skin  friction  and  momen- 
tum thickness  as  parameters.  Additional  parameters  required  in  the  suggested  model 
are  heat  transfer,  energy  thickness,  roughness  height  and  a measure  of  a normal 
entropy  gradient.  Roughness  and  entropy  appear  as  known  boundary  conditions, 
leaving  four  dependent  parameters  to  be  found  from  the  solution  of  the  four  equa- 
tions: velocity  profile,  temperature  profile,  integral  momentum,  and  integral 

energy  equations,  The  advantages  of  this  approach  is  the  accuracy  derived  in 
handling  nonsimilar  boundary  layer  effects  (such  as  step  change  in  surface  mass 
flux,  wall  temperature,  etc.)  relative  to  th-  computationally  simpler  approach  in 
(1)  above.  The  wall  momentum  and  heat  flu.,  variations  are  then  derived  quantities 
similar  to  those  obtained  in  exact  solutions,  but  obtained  more  simply  than  the 
exact  solutions. 
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For  either  of  the  above  approaches,  the  degree  of  success  is  in  part  dependent 
upon  the  availability  of  good  boundary  layer  data  to  be  used  for  correlation  purposes, 
that  is  data  to  generate  (for  example)  wall  shear  laws  as  required  in  (1)  above  or  to 
generate  (for  example)  boundary  layer  wake  laws  as  required  in  ( 2 ) ‘'nve.  These  kinds 
of  data  are  difficult  to  obtain  experimentally,  and  costs  are  usually  prohibitive, 
Accordingly,  exact  codes  such  as  BLIMP  serve  as  an  excellent  resource  for  generating 
the  requisite  data.  However,  the  numerics  must  be  clean  and  the  physics  must  be  repre- 
sented accurately.  The  somewhat  tortuous  variation  of  BUMP  heat  transfer  data  in 
Figure  4-6  is  believed  to  be  due  to  identified  numerical  difficulties.  Since  BLIMP  is 
considered  to  be  an  important  contributor  to  the  success  of  possible  improvements  of 
integral  techniques,  both  the  numerical  problems  and  insufficiencies  of  the  physics 
of  vortical  flow  modeling  that  are  known  to  exist  in  BLIMP  are  discussed  in  the  next 
subsection,  along  with  suggested  means  for  improving  BLIMP. 

4.4  POSSIBLE  EXTENSIONS  OF  BLIMP  TURBULENCE  MODELS  TO  INCLUDE  VORTICAL  EDGE  CONDITIONS 

Previous  sections  have  shown  the  need  for  improvement  of  the  turbulent  forced  convec- 
tion prediction  methodology  to  provide  more  accurate  representations  of  vortical  layer  effects. 

A purely  empirical  approach  has  the  Inherent  difficulty  Of  uncertainty  of  applicability  to 
conditions  which  depart  significantly  from  the  conditions  which  provide  the  empirical  data. 

If  feasable,  it  is  more  desireable  to  incorporate  empirical  information  into  some  form  of  a 
fundamental  theoretical  model,  thereby  minimizing  the  restrictions  on  the  pertinence  of  the 
empirical  information.  The  BLIMP  computer  code  represents  a valid  and  fundamental  theoretical 
modpl  which  is  ideally  suited  for  examining  vortical  layer  effects  in  a very  basic  way.  Its 
end  uses  are:  (1)  the  provision  of  data  pertinent  to  vortical  layer  flows,  and  (2)  the  eval- 

uation of  the  utility  of  correlations  which  are  derived  from  BLIMP  output  and  are  used  in 
simpler  procedures.  Consequently,  extensions  of  the  BLIMP  code  to  Include  vortical  edge  con- 
ditions in  turbulent  flow  are  considered  in  this  subsection,  along  with  a brief  presentation 
of  certain  of  the  elements  currently  existing  in  the  code. 

The  BLIMP  code  solves  to  any  desired  degree  of  accuracy  the  full  boundary  layer  con- 
servation equations  (species,  mass,  momentum,  and  energy)  Including  the  effects  of  kinetically 
controlled  equilibrium  chemistry  (homogeneous  and  heterogeneous),  equal  or  unequal  diffusion 
coefficients,  thermal  diffusion  and  diffusion  thermal  (the  Dufour  effect).  The  solution  pro- 
cedure utilizes  an  efficient  Newton-Raphson  iteration  procedure  to  solve  simultaneously  at 
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nodal  points  across  the  boundary  layer  the  algebraic  equations  formed  by  the  linearized  con- 
servation equations  and  spline  curve  fits.*  (See  Reference  4-13  for  a full  description  of 
this  technique.)  A special  option  has  been  built  into  BLIMP  which  allows  for  the  treatment 
of  edge  normal  velocity  gradients  and  "entropy  swallowing".  The  combined  problem  of  normal 
velocity  gradients  and  "entropy  swallowing"  will  be  subsequently  termed  the  vorticity  problem. 
Given  the  coordinates  of  the  bow  shock  and  the  surface  pressure  distribution,  BLIMP  sets  up 
a table  of  boundary  layer  station  inviscid  properties  and  their  gradients  as  a function  of 
streamfunction.  During  the  boundary  layer  solution  procedure,  the  streamfunction  Is  utilized 
to  couple  the  inviscid  to  the  boundary  layer  flow,  so  that  boundary  layer  edge  properties  and 
gradients  are  matched  to  the  inviscid  values  of  every  iteration.  Although  this  procedure  in- 
cludes the  important  effects  of  edge  normal  gradients  and  "entropy  swallowing",  it  is  not  a 
formally  consistant  approach  since  displacement  effects,  lateral  pressure  gradient  effects 
and  other  higher  order  viscous  effects  are  not  included  (see  for  example,  Reference  4-14  for 
a discussion  of  higher  order  effects  in  boundary  layers).  For  missile  nosetip  conditions,  the 
vortical  layer  is  the  most  significant  of  second  order  effects,  and  the  approach  taken  in 
BLIMP  is  appropriate  for  these  problems.  However,  in  the  case  of  turbulent  flow  with  edge 
vorticity,  problems  arise  in  the  implementation  of  the  currently  used  phenomeiological  tur- 
bulent eddy  viscosity  model  creating  spurious  deviations  from  the  experimental  results.  This 
is  evident  in  Figure  4-6  where  BLIMP  results  using  the  Kendall  turbulence  model  and  tne  vor- 
ticity option  which  conventional  computational  constraints  (velocity  constraint)  are  compared 
with  experimental  data.  The  unnatural  variation  observed  in  this  figure  does  not  appear  for 
turbulent  flows  without  vorticity.  Also,  from  Figure  4-6  it  is  observed  that  neither  the 
(.ebeci/Smith  or  Bushnell/Beckwith  turbulence  models  improve  the  agreement  with  data  over  that 
of  the  Kendull  model.  The  root  of  the  problem  and  a possible  course  for  its  solution  can  be 
found  in  the  examination  of  the  phenomenological  turbulence  models  and  their  implementation 
into  BLIMP. 

All  of  the  above  mentioned  turbulence  models  make  use  of  the  eddy  viscosity  concept 
where  the  turbulent  Reynolds'  stresses  are  written  as: 
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Higher  order  curve  fits  between  nodes  are  utilized  to  allow  matching  of  edge  property  gra- 
dients as  well  as  the  property  values.  This  is  required  to  include  second  order  viscous  velo- 
city gradient  ( 1 . e . , vorticity)  effects  in  a first  order  boundary  layer  treatment. 
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where  t is  the  eddy  viscosity, 
m 

Similarly  for  heat  conduction: 
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where  is  the  eddy  conductivity  and  is  related  to  the  eddy  viscosity  throuqh  the  turbulent 
Prandtl  number  (a  value  of  0.9  is  used  for  turbulent  Prandtl  and  Schmidt  numbers,.  All  of 
the  models  make  use  of  the  Prandtl  mixing  length  hypothesis  in  which  is  is  assumes  that 

c = (4-27) 


where  c,  is  some  mixing  length  and  vt  is  the  turbulent  velocity.  The  differences  between  the 
models  is  in  the  manner  in  which  t and  vf  are  formulated.  In  the  Kendall  and  Cebeci  models 
the  boundary  layer  is  divided  into  an  inner  “wall"  and  and  otter  'wake"  region.  For  the  "wall" 
region  in  both  cases 


v 


t 


(4-28) 


where  is  formulated  so  as  to  vanish  at  the  wall  and  Increase  with  distance  f-om  the  wall. 

In  these  models  the  eddy  viscosity  is  assumed  to  be  a function  of  local  mean  qjantities  (i.e., 
shear  and  distance  from  the  wall)  rather  than  global  quantities  (e.g.,  external  stream  velo- 
city and  boundary  layer  thickness).  This  is  in  keeping  with  the  concept  that  turbulence  in 
the  ''wall"  region  adjusts  rapidly  to  local  conditions  and  does  not  have  much  'historical" 
dependence  on  upstream  conditions.  On  the  other1  hand,  the  outer  "wake"  region  turbulence  has 
historical  dependence  and  therefore  a model  based  on  detailed  local  mean  quantities  doesn't 
give  any  advantage  over  the  more  simplified  global  approach.  Recognizing  this,  as  well  as 
the  fact  that  turbulent  eddy  viscosity  is  nearly  constant  across  the  "wake"  region,  Kendall 
and  Cebeci  applie-'  the  model 
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where 


is  the  velocity  defect  thickness.  The  f LIMP  procedure  incorporates  a coordinate  stretching 
system  which  enables  all  but  the  outermost  nodes  to  be  always  within  the  boundary  layer,  in- 
dependent of  boundary  layer  thickness.  The  coordinate  system  is  forced  to  stretch  with  the 
boundary  layer  by  requiring  one  of  the  nodes  (the  match  node  as  controlled  by  input  designa- 
tion) to  always  have  a velocity  which  is  a constant,  fraction  of  the  edge  velocity.  Thus, 
the  outermost  node  is  required  to  have  identically  the  "edge"  velocity,  and  it  is  located  at 
a dis  ance  from  the  wall  which  scales  with  the  distance  from  the  wall  of  the  1 latch"  node. 
Accordingly,  the  value  of  6 is  not  determined  in  BLIMP.  Thus,  the  value  of  ue  that  BLIMP 
finds  is  the  inviscid  velocity  that  would  exist  at  the  stream  function  of  the  outermost  node. 
Except  for  the  edge  iitermittency,  the  eddy  viscosity  is  assumed  constant  across  the  "wake", 
being  dependent  on  the  global  parameters  of  edge  velocity  and  velocity  defect  thickness. 

The  quantities  and  above  are  not  known  individually,  but  their  product  is  known  from 
experiment  (CjC2  = 0.018)  and  only  their  product  need  be  known  to  implement  the  computations . 

Unlike  the  Kendall  and  Cebeci  models,  the  Bushnell /Beckwith  model  assumes  that  a single 
expression  for  eddy  viscosity  applies  throughout  the  boundary  layer.  In  this  expresison 
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where  H is  a function  of  distance  from  the  wall,  boundary  layer  thickness,  and  interim ttency 
(as  a function  of  y/A)  such  that  mixing  length  and  eddy  viscosity  decay  to  zero  at  the  boun- 
dary layer  edge. 

For  irrotational  external  flows  (no  vorticity),  all  of  the  models  are  in  reasonable 
agreement  amongst  themselves  as  well  as  with  the  data  (c.f.  Reference  4-15).  For  rotational 
external  flows,  it  can  be  concluded  from  the  discussion  above  of  the  edge  condition  in  BLIMP 
that  evaluations  ->f  Equations  (4-29  through  (4-31)  are  all  sensitive  to  the  location  of  the 
boundary  layer  "edge".  From  these  considerations  evolved  the  approach  of  designating  enthalpy 
of  the  "match"  node  rather  than  velocity,  since  at  least  the  edge  stagnation  enthalpy  is  not 
a function  of  edge  location,  unlike  the  "edge"  velocity. 


Thus,  regardless  of  the  physical  applicability  of  the  Cebeci/Smith  or  Kendall  eddy 
viscosity  models  to  flows  with  entropy  swallowing,  it  is  evident  that  the  current  implementa- 
tion of  these  models  restricts  their  validity  to  flows  with  no  external  vorticity.  In  the  fol- 
lowing paragraphs,  the  physical  applicability  of  the  eddy  viscosity  model  to  flows  with  vor- 
ticity is  briefly  discussed.  This  is  followed  by  a discussion  of  the  means  of  correctly  im- 
plementing the  model  for  flows  with  vorticity. 

Before  discussing  turbulent  transport  models  it  should  be  mentioned  that  the  boundary 
layer  concept  is  only  valid  (and  only  approximately  so)  up  to  a luvel  of  free  stream  vorticity 
which  is  a fraction  of  that  in  the  boundary  layer  itself  (i.e.,  "free  stream"  shear  stresses 
that  exist  due  to  finite  viscosity  should  be  small).  If  the  vorticity  is  larger  than  the 
boundary  layer  value,  no  matter  what  turbulent  transport  model  is  utilized,  boundary  layer 
predictions  of  ht.at  flux  will  not  be  adequate,  Therefore,  the  discussion  is  restricted  to 
turbulence  models  which  have  been  applied  to  flows  with  no  edge  vorticity  and  which  can  be 
easily  extended  to  the  cases  with  moderate  or  small  external  stream  vorticity. 

As  previously  mentioned,  the  turbulence  near  the  surface  (i.e.,  "wall"  region)  adjusts 
rapidly  to  local  conditions  and  is  insensitive  to  upstream  and  boundary  layer  edge  events. 

Since  the  external  vorticity's  impact  on  the  boundary  layer  is  primarily  in  altering  edge 
conditions  and  creating  a historical  dependence,  currently  used  "wall"  turbulent  models  should 
be  adequate.  For  the  “wake"  region  of  the  boundary  layer,  external  vorticity  will  probably 
modify  the  classical  generation,  convection,  diffusion,  and  dissipation  of  turbulence.  Pre- 
sently, there  are  no  experimenta1  or  theoretical  results  which  indicate  the  manner  in  which 
these  processes  are  affected  by  vorticity.  However,  it  may  be  conjectured  that  the  external 
vorticity  would  enhance  the  production  (i.e.,  du/dy  t 0 at  edge)  of  large  scale  turbulent 
motions.  Since  the  impact  is  through  the  edge,  the  classically  observed  edge  intermittency 
would  be  modified  as  a result  of  stream  vorticity.  A detailed  accounting  for  these  effects 
might  be  possible  using  a two  equation*  detailed  model  of  turbulence.  However,  the  added  dif- 
ficulty in  solving  these  equations  and  the  lack  of  knowledge  concerning  the  empirical  coeffi- 
cients would  not  justify  this  approach  at  the  present  time.  A more  straight  forward  and 
simpler  approach  Is  to  assume  that  the  external  vorticity  has  a global  effect  on  the  "wake" 
portion  of  the  boundary  layer  and  therefore  existing  phenomenaloglcal  models  can  be  utilized 
if  the  empirical  coefficients  and  means  of  implementation  are  adjusted  accordingly.  As 
* 

Differential  equations  describing  turbulent  kinetic  energy  and  length  scale  would  be  added 

to  the  basic  conservation  equations. 


previously  mentioned  no  advantage  is  gained  In  an  approach  which  depends  on  local  quantities 
in  the  outer  portion  of  the  boundary  layer.  Therefore,  a "wall"  and  “wake"  law  phenomenolog- 
ical approach  similar  to  those  of  Kenda'l  and  Cebeci  can  be  adapted  for  use  in  the  case  with 
external  vorticity.  As  in  the  case  w' chout  vorticity,  an  expression  for  vt  and  l must  be 
found  to  model  the  eddy  viscosity  in  the  wake  portion  of  the  boundary  layer,  and  these  ex- 
pressions should  degenerate  to  their  conventional  forms  in  the  absence  of  a vortical  layer. 

The  most  straight-forward  modification  to  achieve  this  end  is  to  retain  Equations  (4-29) 
and  (4-30)  as  they  stand,  but  modify  the  velocity  defect  thickness  to  the  following  form: 

6i'/=i  ['  ■Jr]"'  (4-331 

where  u.  is  the  inviscid  velocity  that  would  exist  at  the  local  boundary  layer  stream  function. 
(The  actual  value  of  velocity  at  the  local  stream  function  is  u.)  Note  that  u/u.  = 1 every- 
where outside  of  the  boundary  layer,  independent  of  the  magnitude  or  variation  of  u^,  so  that 
is  insensitive  to  the  upper  limit  on  the  integration,  providing  orly  that  u is  approxinately 
equal  to  u-  at  the  upper  limit.  Note  that  for  conventional  flows,  i...  = ug  so  that  Equation 
(4-31)  is  a restricted  case  of  the  more  general  expression,  Equation  (4-33). 

Obviously,  the  evaluation  of  Equation  (4-33)  is  still  sensitive  to  the  determination 
of  ufi,  the  boundary  layer  edge  velocity.  Since  there  is  ambiguity  concerning  the  location  of 
5,  ufi  will  be  sensitive  to  this  choice  when  du/dy  / 0.  However,  it  can  be  noted  that  the  eddy 
viscosity  is  completely  insensitive  to  this  choice  when  Equations  (4-27),  (4-29),  (4-30),  and 
(4-33)  are  combined. 


6 

= lvt  = ClC2ue6i  = C1C2  f ^'i  ” ^ dy  (4-34) 

o 

In  Equation  (4-34),  the  precise  value  of  6 is  unimportant,  since  u * u^  regardless  of 
location  in  the  outer  reaches  of  the  boundary  layer.  Thus,  the  proposed  modification  at  lease 
obviates  certain  numerical  difficulties  that  exist  with  the  present  procedure,  which  is  sen- 
sitive to  the  value  of  6.  This  is  particularly  true  for  cases  where  u^  does  not  vary  mono- 
tonically  with  y,  as  happens  when  the  boundary  layer  edge  streamline  passes  through  an  inflec- 
tion in  the  shock  (the  simple  biconic  geometry  of  Figure  4-6  has  such  an  inflection  - the  re- 
duction in  computed  heat  flux  at  S/R  0.75  using  the  velocity  constraint  coincides  with  an 


4-37 


abrupt  change  In  sign  of  the  normal  entropy  gradient  near  the  boundary  layer  edge  caused  by 
this  shock  inflection). 

Besides  the  proposed  alternation  of  the  "wake"  law,  it  might  also  be  important  to  al- 
ter the  edge  intermittency  model.  Recognizing  that  free  stream  vorticity  might  enhance  pro- 
duction of  turbulence  in  the  outer  portion  of  the  boundary  layer,  a slower  eddy  viscosity 
decay  to  possibly  even  a nonzero  value  might  be  a reasonable  condition  to  apply  for  flows 
with  vorticity. 

In  summary,  it  is  believed  that  existing  phenomenological  eddy  viscosity  models  can 
be  suitably  modified  to  apply  to  flows  with  vorticity.  The  validity,  flexibility  and  accu- 
racy of  such  models  can  only  be  established  through  comparisons  with  experimental  data.  In 
this  regard,  the  need  for  detailed  boundary  layer  survey  data  for  these  types  of  flows  is 
apparent. 

4.5  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 

The  requi renents  of  the  computational  methodology  for  turbulent  convective  transfer 
phenomena  for  use  in  shape  change  codes  have  been  described  herein.  From  this  discussion, 
the  desirability  of  the  use  of  integral  boundary  layer  techniques  is  suggested,  based  primarily 
on  computation  speed  requirements. 

The  integral  technique  for  the  turbulent,  nonblown,  roughwall  boundary  layer  evaluation 
that  is  used  in  the  EROS  shape  change  code  is  presented.  In  addition,  certain  of  the  details 
of  the  more  refined  exact  turbulent  boundary  solution  methodology  employed  in  the  BLIMP  code 
are  presented.  Results  from  bo,.,.  of  these  techniques  have  been  compared  with  appropriate 
experimental  data.  Certain  techniques  were  tried  to  improve  the  comparisons  of  predictions 
with  data,  and  cc-'tain  desirable  modifications  to  the  procedures  have  been  identified. 

Efforts  'o  improve  results  by  considering  variable  Reynolds1  analogy  factors  in  the 
EROS  code  as  derived  from  exact  BLIMP  calculations  were  generally  successful.  Also,  attempts 
to  improve  predictions  through  an  averaging  scheme  for  the  swallowed  entropy  layer  were  appar- 
ently successful  relative  to  ballistic  range  data  at  high  Mach  number,  but  not  successful  ior 
the  M = 5 wind  tunnel  data.  The  basic  EROS  scheme  without  modifications  predicts  the  data  in 
the  high  heat  flux  regions  of  the  bodies  considered  within  about  i 25  percent  for  the  M = 5 
experiments.  The  basic  inegral  scheme  can  be  inaccurate  for  high  Mach  number  environments 
when  a strong  entropy  swallowing  effect  occurs  on  the  forecone.  In  general  the  prediction 
lies  above  the  data.  In  some  cases  considered  the  local  heat  flux  is  as  much  as  a factor  of 
t.wo  high. 
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The  basic  EROS  turbulent  boundary  layer  technique  consists  of  solution  of  a simplified 
form  of  the  momentum  integral  equation.  The  similarity  between  the  momentum  integral  equation 
and  the  energy  integral  equation  is  noted.  It  is  shown  that  the  EROS  procedure  does  not  ac- 
count for  the  “acceleration"  terms  in  either  of  these  equations.  However,  it  is  shown  that 
the  present  procedure  is  in  fact  a solution  to  the  energy  integral  equation  for  an  isothermal 
wall,  as  a consequence  of  the  simplification  that  H = -1  in  the  integral  momentum  equation. 

The  poor  agreement  of  the  basic  EROS  scheme  with  the  ballistic  range  data  is  attributed 
to  both  poor  representations  of  boundary  layer  shape  factors  in  the  boundary  layer  entropy 
swallowing  relation  and  excessively  large  estimates  of  boundary  layer  momentum  thickness  which 
are,  in  effect,  assumed  to  be  equal  to  boundary  layer  energy  thicknesses.  In  addition,  com- 
puted momentum  thicknesses  are  too  large  as  a consequence  of  the  use  of  the  composite  heating 
relationship.  It  is  also  expected  that  entropy  layer  effects  will  be  poorly  represented  in 
flight  cases  where  the  transpiration  effects  which  are  not  presently  accounted  for  influence 
both  the  boundary  layer  parameters  and  the  entropy  swallowing  relation  (although  these  ef- 
fects are  somewhat  compensatory). 

The  lack  of  improvement  of  comparisons  of  EROS  predictions  with  data  using  variable 
Reynolds'  analogy  factors  is  ascribed  to  the  postulate  that  It  is  the  integral  energy  boundary 
layer  equation  that  is  really  being  solved.  As  such,  the  solution  has  no  need  for  inclusion 
of  the  variable  Reynolds'  analogy  factor,  except  to  perhaps  estimate  by  the  inverse  of  the 
Reynolds'  analogy  factor  the  friction  factor  from  the  derived  Stanton  number.  It  is  suggested 
that  the  Reynolds'  analogy  factor  that  should  be  employed  is  that  derived  from  BLIMP  calcula- 
tions for  a flat  plate  (Rj  = Pr'V3). 

Of  the  second  order  boundary  layer  effects  that  might  be  considered  for  shape  change 
codes,  it  is  believed  that  the  entropy  layer  effect  is  most  important,  particularly  for  the 
high  flight  Mach  numbers  where  the  entropy  layers  are  swallowed  in  the  forward  regions  of  the 
body.  The  BL’Mr  code  exhibits  some  computation  difficulties  with  vortical  layers  when  attempts 
are  made  to  apply  conventional  turbu1ent  boundary  layer  wake  and  intermi ttency  models  to 
these  unconventional  cases.  The  BLIMP  code  was  exercised  for  M = 5 wind  tunnel  data  employing 
three  different  wake  formulations,  including  the  methods  of  Kendall,  Cebeci/Smith,  and  Bushnell/ 
Beckwith.  All  models  predicted  the  shape  of  the  heat  transfer  distributions,  and  the  Kendall 
model  predicts  the  level  as  well.  The  Kendall  model  also  does  reasonably  well  at  high  Mach 
number,  and  the  means  to  improve  the  BUMP  methodology  in  order  to  obviate  the  need  for  arti- 
ficial constraints  for  all  applicable  vortical  flows  have  been  identified. 


The  need  to  extend  the  integral  method  to  solution  of  both  the  momentum  and  energy 
equations  has  been  suggested.  The  momentum  boundary  layer  solution  is  needed  to  characterize 
transitional  phenomena  and  to  provide  accurate  mass  balances  in  the  boundary  layer.  iMs 
solution  should  then  yield  accurate  boundary  layer  state  and  edge  conditions  for  evaluation 
of  the  energy  boundary  layer,  thereby  improving  the  turbulent  convective  heat  and  mass  flux 
(by  analogy)  predictions. 

Two  levels  of  improvement  to  the  integral  boundary  layer  solution  methodology  are  con- 
sidered. The  simpler  method  from  a computation  time  standpoint  requires  correlation  of  wall 
shear  and  heat  flux  data  to  account  for  wall  roughness,  transition,  transpiration,  free 
stream  acceleration,  wall  temperature  gradients,  and  vortical  layer  effects.  The  more  so- 
phisticated method  requires  correlation  of  boundary  layer  velocity  and  stagnation  enthalpy 
distribution  data  in  terms  of  a “law  of  the  wall11  and  “law  of  the  wake  to  account  for  the 
same  effects  considered  in  the  simpler  approach  above.  The  advantages  to  this  latter  approach 
is  the  expected  increased  accuracy  obtained  without  the  need  for  a-priori  specifications  of 
how  the  wall  flux  laws  vary  with  the  phenomena  that  need  to  be  considered  (roughness,  trans- 
piration, etc.).  The  expected  disadvantage  is  -'ncreased  computation  time. 

The  success  of  either  of  the  integral  methods  above  are  somewhat  dependent  upon  the 
availability  of  good  boundary  layer  data  to  be  used  for  correlation  purposes.  The  BLIMP  code 
with  changes  as  recommended  herein  should  serve  as  an  excellent  resouce  for  generation  of 
these  data,  since  the  code  capabilities  are  quite  general.  Experimental  boundary  layer  data 
with  vortical  boundary  layer  edge  conditions  are  currently  nonexistent,  but  would  be  very 
useful  in  terms  of  verifying  or  providing  insight  for  improvement  of  the  BLIMP  models.  Ef- 
forts to  plan  and  conduct  these  experiments  are  recommended. 
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SECTION  5 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  primary  conclusions  and  recommendations  resulting  from  the  invl'cld  flow,  smooth 
wall  nosetip  transition,  and  turbulent  heat  transfer  tasks  are: 

• Inviscid  Flow 

— The  thin  shock  layer  approach  was  demonstrated  to  yield  good  potential  for  pre- 
dicting accurate  shock  shapes  when  the  surface  pressure  distribution  is  spec- 
ified. Additional  dcvelonment  in  the  solution  procedure  is  necessary,  and  is 
recommended,  to  obtain  a procedure  which  is  adequate  for  inclusion  in  the  shape 
change  codes. 

• Smooth  Wall  Nosetip  Transition 

- On  smooth  blunt  nosetips,  the  boundary  layer  ma^  remain  laminar  to  Quite  high 
Reynolds  numbers  (i.e.,  comparatively  low  altitudes). 

— Transition  on  smooth  surfaces  is  sensitive  to  mihor  surface  irregularities 
and/or  variations  in  atmospheric  conditions. 

— Transition  on  smooth  reentry  nosetips  will  be  difficult  to  predict  accurately 
and  may  be  fundamentally  uncertain  because  of  atmospheric  variations. 

Therefore,  if  smooth  materials  are  important  in  fulfilling  mission  requirements, 
studies  should  be  conducted  to  determine  the  sensitivity  of  nosetip  performance 
to  uncertainties  in  transition  altitude. 

• Turbulent  Heat  Transfer 

- Theoretical  considerations  indicate  that  the  thickness  which  is  obtained  by 
solving  the  simplified  momentum  equation  is  probably  more  characteristic  of  the 
energy  than  the  momentum  layer  thickness.  Preliminary  indications  are  that  use 
of  this  thickness  in  the  entropy  layer  effects  mass  balance  can  result  in  over- 
prediction  of  turbulent  heat  and  mass  transfer  rates. 
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In  the  transitional  regime,  computed  momentum  thicknesses  are  believed  to  be 
too  large  as  a result  of  the  use  of  the  composite  heating  relationship. 

The  reconmended  turbulent  Reynolds  analogy  factor  for  use  with  the  current 
procedure  is  Pr  n . 

Although  the  present  technique  fails  to  predict  the  ballistic  range  data  con- 
sidered herein,  adequate  predictions  for  flight  conditions  can  be  realized  by 
extending  the  integral  solution  to  include  both  the  momentum  and  energy  equa- 
tions. The  momentum  boundary  layer  solution  is  necessary  to  characterize  tran- 
sition and  to  provide  accurate  mass  balances  in  accounting  for  entropy  layer 
effects.  The  energy  boundary  layer  solution  is  necessary  to  adequately  account 
for  convective  heat  and  mass  transfer  effects. 

Experimental  and  numerical  data  are  needed  to  adequately  extend  and  validate  the 
heat  and  mass  transfer  prediction  procedures. 
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